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Abstract 

■V 

The  a  priori  hypothesis  that  soli -lunar  atmospheric  tides  can 

~~  V> 

Induce  an  element  of  the  variability  in  mean  monthly  climatic  data  is 
tested  using  precipitation  (1895-1960),  station  pressure  (1895-1960) 
and  300  mb  height  data  (1950-1977)  for  India  during  the  monsoon  onset 
month  of  June.  Since  the  monsoon  onset  occurs,  on  average,  very  late 
In  the  month  In  northern  India,  this  region  Is  very  sensitive  to 
small  departures  from  normal  In  the  onset  date.  These  departures  can 
be  seen  In  the  year  to  year  fluctuations  of  June  rainfall  amounts. 

To  test  the  hypothesis  involving  lunar  atmospheric  tides,  a  time 
series  of  mean  monthly  lunar  tidal  potential  for  June  (1895-1960)  was 
calculated."  A  perlodogram  of  this  time  series  is  dominated  by  two 
frequencies,  .0537  year"*  and  .263  year"*,  which  explain  almost 
100  percent  of  the  variance.  The  frequency  spectrum  of  the  first 
eigenvector  of  June  precipitation  Is  dominated  by  these  two 
frequencies  with  the  frequency  .0537  year"*  being  shifted  somewhat 
lower.  This  frequency  shift  was  caused  by  a  nonlinearity  In  the 
atmospheric  response  to  the  forcing.  It  was  found  that  the  amplitude 
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of  the  response  was  not  constant  in  time  but  dependent  upon  the 
temperature  of  the  Northern  Hemisphere  which  determines,  to  a  large 
extent,  the  intensity  of  the  monsoon  circulation.  When  the 
circulation  Is  strong  (l.e.,  when  the  Northern  Hemisphere  is  warm), 
the  effects  of  lunar  atmospheric  tides  on  climatic  variability  are 
minimal  and  vice  versa.  Thus,  the  atmosphere  is  more  susceptible  to 
lunar  forcing  at  certain  times  than  at  others  making  the  atmospheric 
response  nonlinear.  The  ramifications  of  this  nonlinear  response  are 
discussed  in  depth. 

i'r  0  > The  results  of  this  research  showed  that  there  is  an  important 
periodic  component  (on  the  order  of  one  standard  deviation  during 
hemispherical ly  cold  periods)  in  the  interannual  variability  of 
pressure,  precipitation,  and  upper  air  winds  for  the  Indian  monsoon 
onset  month  of  June  at  the  two  dominant  lunar  frequencies  and  that 
the  atmospheric  response  is  In  phase  with  the  forcing.  This  strongly 
suggests  that  lunar  atmospheric  tides  do  produce  an  element  of  the 
observed  Interannual  variability* 
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CHAPTER  1 


Introduction 


1.1  Background. 

The  planetary  monsoon  circulation  Is  a  dominant  element  of  the 
large-scale  atmospheric  response  to  the  annual  solar  cycle.  As  the 
Northern  Hemisphere  warms  In  the  spring,  the  expanse  of  the 
westerlies  contracts,  the  subtropical  highs  move  northward  and  the 
onset  of  the  monsoon  begins.  Since  this  planetary  phenomenon,  which 
Includes  the  southwest  summer  monsoon  In  India,  results  from  thermal 
Inhomogeneities  at  the  earth's  surface  created  by  the  distribution 
of  the  continental  land  mass  and  oceans,  the  large-scale  features  of 
the  monsoon  circulation  are  similar  from  one  year  to  the  next. 
However,  significant  Interannual  variability  of  rainfall  exists  In 
the  monsoon  region.  For  example,  very  heavy  rains  fell  In  India 
during  the  1917  monsoon  whereas  the  1918  monsoon  season  produced  one 
of  the  worst  droughts  of  the  last  century.  A  reconstructed  history 
of  monsoon  precipitation  (Bryson  and  Swain,  1981)  shows  that  the 
monsoon  has  not  been  a  constant  and  reliable  phenomenon,  and  this 
variability  will  almost  certainly  continue  In  the  future.  With  the 
econonty  of  India  being  so  fragile  that  one  drought  Is  a  sufficient 
factor  alone  to  disrupt  the  country.  It  Is  Imperative  that  we 
develop  the  capability  to  forecast  the  year  to  year  fluctuations  In 
precipitation  If  anyone  Is  to  plan  rationally  for  the  future  of  the 
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June  Is  the  month  of  the  monsoon  onset  over  most  of  the  Indian 
subcontinent.  The  onset  can  be  viewed  as  a  front  (Das,  1968)  which, 
on  the  average,  enters  the  southeastern  tip  of  India  around  1  June 
and  reaches  the  Pakistani -Indian  border  by  1  July  (Fig.  1.1).  Prior 
to  monsoon  onset,  India  Is  basically  dry  and  hot  except  for 
scattered  pre-monsoon  showers.  When  the  monsoon  front  arrives  there 
Is  a  tremendous  burst  of  precipitation.  However,  once  the  monsoon 
front  Is  past,  the  cloud  tops  are  generally  not  very  high  and 
drizzle  to  light  rain  Is  the  predominant  form  of  precipitation  with 
Intervals  of  heavier  rainfall  and  deeper  convection  associated  with 
monsoon  depressions  from  the  Bay  of  Bengal. 

There  Is  considerable  variability  In  the  date  of  onset  from  year 
to  year.  For  example,  along  the  west  coast  south  of  20°N  the 
standard  deviation  of  the  onset  Is  about  seven  days  (Rao,  1976).  It 
Is  the  variability  In  onset  date  which  greatly  affects  the  June 
precipitation  pattern  In  northern  India,  where  the  advance  of  the 
monsoon  front  begins  to  slow  down.  In  this  region  small 
displacements  of  the  monsoon  front  with  Its  associated  large 
rainfall  gradient  will  markedly  affect  the  amount  of  rainfall 
received  In  June.  If  the  monsoon  front  Is  late,  a  significant 
percentage  of  the  opportunity  for  rain  In  the  calendar  month  of  June 
Is  lost,  because  the  pre-monsoon  days  are  drier  than  the  normal 
post-monsoon  "burst"  days.  Thus,  small  anomalies  In  the  monsoon 
front's  normal  position  can  be  seen  In  the  year  to  year  fluctuations 
of  June  rainfall  amounts  at  stations  In  northern  India.  For 
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example,  If  the  normal  date  of  arrival  of  the  monsoon  is  25  June,  a 
two  day  delay  In  arrival  time  represents  a  40  percent  reduction  In 
the  duration  of  June  monsoon  rains.  By  contrast,  at  stations  In 
southern  India,  where  the  monsoon  normally  arrives  In  Hay  or  early 
June,  a  two  day  delay  In  arrival  of  the  rains  would  have  a  small 
percentage  effect  on  June  rainfall.  Furthermore,  In  northern  India 
the  monsoon  front  moves  slowly,  on  the  order  of  25-50  km  per  day  In 
the  mean,  so  that  quite  small  spatial  departures  from  the  normal 
position  repiesent  significant  advances  or  delays  In  arrival  time. 

Once  the  monsoon  Is  established  over  India,  rainfall  Is  the 
result  of  a  very  complex  pattern.  There  are  two  related  systems. 

One  brings  rains  to  the  west  coast  with  the  strong  southwesterly 
onshore  flow,  and  the  other  brings  rains  to  the  Gangetlc  Plain  where 
most  of  the  rain  In  July  and  August  Is  associated  with  monsoon  lows 
and  depressions  which  develop  In  the  northern  Bay  of  Bengal  and 
travel  westward  across  northern  India.  On  top  of  these  large  scale 
features  are  then  superimposed  substantial  effects  due  to  topography 
along  the  western  Ghats  and  the  Himalayan  foothills  which  tend  to 
obscure  the  Influence  of  the  large-scale  circulation  features. 

Added  to  this  already  complex  pattern  Is  the  movement  of  the  monsoon 
trough.  As  Raghavan  (1973)  points  out,  the  track  of  a  monsoon 
low/depression  and  the  position  of  the  monsoon  trough  are  related  - 
with  a  more  southerly  position  being  favorable  for  rain  over  most  of 
the  country  and  a  more  northerly  position  against  the  Himalayas 
being  associated  with  breaks  In  the  monsoon. 
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The  Indian  sunner  monsoon  has  been  the  subject  of  extensive 
study  for  many  years.  For  example,  Krishnamurtl  has  studied  the 
elements  of  the  monsoon  circulation  In  detail  (e.g.,  Krishnamurtl, 
1971;  Krishnamurtl  et  al.,  1973;  Krishnamurtl  and  Bhalme,  1976). 
Others  have  performed  observational  studies  on  the  Interannual 
variability  of  Indian  sunner  monsoon  rainfall.  Hahn  and  Shukla 
(1976)  showed  an  apparent  association  with  Eurasian  snow  cover  -  the 
more  the  snow  the  less  the  rain.  Shukla  and  Mlsra  (1977)  and  Weare 
(1979),  among  others,  have  studied  the  effects  of  sea  surface 
temperatures  on  the  amount  of  Indian  rainfall  with  differing 
results.  Still  others  have  used  numerical  models  to  try  and 
understand  the  monsoon.  For  example,  Hahn  and  Manabe  (1975),  using 
a  general  circulation  model  showed  the  Importance  of  the  Himalayan 
mountains  In  helping  to  extend  the  monsoon  climate  farther  north 
onto  the  Aslan  continent.  Shukla  (1975)  modelled  the  effect  of 
Arabian  sea  surface  temperatures  on  the  monsoon.  Other  numerical 
studies  have  Included  a  simulation  of  the  Indian  suaner  monsoon  by 
Godbole  (1973)  and  a  model  of  the  seasonal  structure  of  a  simple 
monsoon  system  by  Webster  and  Chou  (1980).  Matsuno  (1966)  and  Gill 
(1980)  have  studied  the  atmospheric  response  to  both  symmetric  and 
antisymmetric  heating  about  the  equator  where  Gill's  antisymmetric 
heating  produced  some  of  the  basic  features  of  the  Indian  monsoon 
circulation. 

The  relationship  of  the  Southern  Oscillation  to  droughts  and 
floods  In  India  has  been  the  subject  of  Investigation  ever  since  Sir 


Gilbert  Walker  (1924)  found  a  tendency  for  Increased  precipitation 
In  India  when  the  pressure  Increases  In  the  Pacific  and  lowers  In 
the  Indonesian  region.  Most  recently  Bhalme  et.  al.  (1983)  studied 
the  relationship  between  an  Indian  drought/f 1  ood  area  Index  and  a 
Southern  Oscillation  Index  based  on  April  pressures  and  found 
significant  correlations  between  the  two  Indices.  Using  a  coherence 
spectrum  analysis*  they  concluded  that  the  Southern  Oscillation  Is 
one  possible  mechanism  for  Introducing  a  3-6  year  oscillation  In  the 
drought  area  Index.  The  research  presented  In  the  following 
chapters  will  add  another  mechanism  to  the  possible  causes  of  the 
observed  Interannual  variability  In  the  monsoon  circulation. 

1.2  Research  perspective. 

Climate  can  be  viewed  as  the  state  of  the  earth-atmosphere-ocean 
system  resulting  from  a  set  of  boundary  conditions  Imposed  on  this 
system.  This  climatic  state  then  establishes  an  array  of  weather 
patterns;  and  If  the  climatic  state  changes,  a  different  array  of 
weather  patterns  results.  Thus,  In  order  to  determine  a  climatic 
state  with  Its  attendent  synoptic  patterns  one  must  be  able  to 
specify  the  boundary  conditions,  both  Internal  and  external,  that 
are  Imposed  on  the  atmosphere.  These  are  the  ultimate  causes  of 
climate  and  Include  such  conditions  as  the  energy  output  of  the  sun, 
the  earth's  orbital  characteristics  (e.g.,  time  of  year  of 
perihelion),  pole  and  gravitational  tidal  forcing,  anthropogenic 
heat  and  particulate  output,  and  volcanic  activity.  These  ultimate 
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boundary  conditions  result  In  what  one  might  call  proximal  causes  of 
climate  like  the  amount  of  solar  radiation  reaching  the  earth's 
surface  which  Is  dependent  on  several  of  the  ultimate  causes  given 
above.  If  the  ultimate  causes  vary,  then  the  climate  must  vary 
also.  Thus,  In  order  to  understand  climatic  variability  and  be  able 
to  predict  a  future  climatic  state  we  must  be  able  to  specify  the 
appropriate  boundary  conditions  and  understand  how  they  vary*  Over 
the  last  100  years  some  of  the  boundary  conditions  can  be  assumed 
constant  (e.g.,  earth's  orbital  parameters)  while  others  cannot 
(e.g.,  volcanic  activity,  pole  and  gravitational  tidal  forcing). 
Figure  1.2  shows  a  schematic  which  puts  this  research  effort  Into 
perspective. 

Interannual  variability  of  a  parameter  like  rainfall  amount  can 
be  separated  Into  at  least  two  parts  -  a  periodic  and  a  secular 
component.  The  secular,  or  trend,  component  might  be  associated 
with  the  thermodynamic  state  of  the  atmosphere  and  also  could  be  a 
periodic  component  but  one  with  a  period  near  the  length  of  the  data 
record  or  longer.  A  third  possible  component  might  be  a  stochastic 
one  (l.e.,  the  part  of  climatic  variability  that  Is  currently 
unknown  or  that  Is  truly  random). 

It  Is  Important  to  emphasize  that  there  are  many  causes  of 
climatic  variability  and  change.  The  scope  of  this  research  effort 
has  been  narrowed  to  an  analysis  of  one  of  the  forcing  mechanisms 
that  night  produce  a  periodic  component  In  the  Interannual 
variability  found  In  the  monsoon  onset  month  of  June.  The  research 
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ULTIMATE  CAUSES 

I 

PROXIMAL  CAUSES 
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I 

RESULTING  ARRAY  OF  WEATHER  PATTERNS 
(Monsoon  Circulation) 

1 

INDIAN  SUMMER  MONSOON  RAINFALL 


Figure  1.2  Schematic  of  research  perspective 
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was  conducted  with  the  a  priori  hypothesis  that  soli -lunar 
atmospheric  tides  can  Induce  an  element  of  the  periodic  variability 
In  mean  monthly  climatic  data.  Lunar  tidal  forcing  of  climate  Is 
not  a  new  Idea;  for  example,  Bryson  (1948)  found  that  the 
latitudinal  position  of  the  eastern  Pacific  High  was  highly 
correlated  with  the  equatorward  component  of  the  bi fortnightly  lunar 
tidal  force.  Thirty-five  years  later  there  Is  still  much  research 
needed  In  this  area. 

Chapter  2  will  examine  the  hypothesized  lunar  forcing  by 
developing  a  lunar  tidal  potential  time  series  averaged  over  the 
month  of  June  for  the  years  1895  through  1960.  The  periodic 
components  of  this  time  series  will  be  determined,  and  the  Moon's 
orbital  characteristics  that  produce  this  part  of  the  lunar 
potential  Interannual  variability  will  be  discussed. 

Chapter  3  will  then  test  the  hypothesis  through  a  diagnostic 
study  which  examines  the  response  of  surface  pressure, 
precipitation,  and  300  mb  data  at  the  two  dominant  lunar  forcing 
frequencies.  Campbell  et.  al.  (1983)  found  that  these  two 
frequencies  were  most  Important  In  explaining  the  Interannual 
variability  of  precipitation  In  northern  India.  A  unique  way  of 
defining  how  the  atmosphere  responds  to  the  hypothesized  forcing 
will  be  described  In  the  development  of  a  nonlinear  model  combining 
the  trend  and  periodic  components  of  variability. 

Chapter  4  will  present  the  major  conclusions  of  this  research. 
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CHAPTER  2 

Periodic  Component:  Hypothesized  Forcing 

2.1  Introduction 

Mar\y  studies  of  Indian  summer  monsoon  rainfall  have  found  a 
periodic  component  In  the  data.  Koteswaram  and  Alvi  (1969)  and 
Bhargava  and  Bansal  (1969)  found  quasi-biennial  variations  (QBVs)  In 
both  monsoon  and  annual  rainfall  at  some  stations  on  the  west  coast 
of  southern  India.  Jagannathan  and  Parthasarathy  (1973)  found  a  QBV 
in  the  annual  precipitation  of  other  stations  throughout  India. 
Parthasarathy  and  Mooley  (1978)  in  their  study  of  a  single 
homogeneous  time  series  of  Indian  summer  monsoon  rainfall,  derived 
from  an  areal  weighting  of  over  3000  rain  gauges,  found  a 
significant  cycle  (at  the  5  percent  level)  In  the  QBV  range.  Bhalme 
(1972)  found  the  frequency  of  southwest  monsoon  cyclonic 
disturbances  (1891-1970)  had  a  significant  QBV  of  3.3  to  3.5  years. 
Bhalme  and  Mooley  (1980)  found  a  QBV  frequency  in  the  occurrence  of 
droughts  In  India  and  an  Interval  of  about  20  years  between 
large-scale  floods.  Jagannathan  and  Parthasarathy  (1973)  found 
periodicities  of  approximately  11  years  which  they  related  to  the 
solar  cycle.  Interestingly  enough,  Parthasarathy  and  Mooley  (1978) 
found  no  significant  relationship  between  the  sunspot  cycle  and 
monsoon  rainfall  In  India.  However,  Bhalme  and.  Mooley  (1981)  found 
a  significant  signal  In  their  Flood  Area  Index  (FAI)  at  22  years 


which  they  attributed  to  the  double  (Kale)  sunspot  cycle.  They  used 
a  Fast  Fourier  Transform  analysis  on  88  years  of  data  which 
automatically  puts  band  4  at  a  22  year  periodicity;  whereas,  a 
perlodogram,  which  looks  at  the  spectral  behavior  between  the  bands, 
shows  that  the  spectral  peak  was  actually  around  .048  year'* 

( 1  .e . ,  a  20.8  year  period). 

One  must  ask  what  ultimate  cause  can  produce  periodicities  in 
measured  atmospheric  variables  like  wind,  pressure  and 
precipitation;  for  until  we  can  answer  this  question,  accurate 
climatic  forecasts  of  one  to  five  years  In  advance  will  not  be 
possible.  Bryson  and  Goodman  (1980)  have  shown  that  the  variations 
In  the  amount  of  solar  radiation  reaching  the  ground  have  a 
possible,  non-sunspot  periodic  component.  Feedback  loops  within  the 
earth-atmosphere-ocean  system  may  induce  a  periodicity  In  rainfall 
amounts.  However,  It  seems  hard  to  envision  a  closed  feedback 
system  that  would  oscillate  for  many  years,  unless  It  Is  forced, 
because  of  the  high  correlations  that  would  be  required  between  the 
meteorological  variables.  Another  possible  candidate  Is  the 
periodic,  mechanical  forcing  produced  by  both  the  long  period 
soil-lunar  tides  and  the  "pole  tide"  associated  with  the  Chandler 
wobble.  Studies  by  Bryson  and  Starr  (1977),  Maksimov  (1958)  and 
Maksimov  et  al.  (1967)  have  provided  insight  Into  the  effects  of  the 
pole  tide  on  the  atmosphere. 

Soil-lunar  tidal  forcing  Is  of  particular  interest  since  Its 
periodicities  are  so  highly  predictable.  There  have  been  many 


studies  which  have  related  lunar  periods  to  particular 
meteorological  phenomena  (e.g.,  the  synodic  period1  to  rainfall). 
Brier  and  Bradley  (1964)  examined  16,057  maximum  monthly  24-hour 
precipitation  occurrences  in  the  U.S.  from  1900-1949.  By  relating 
the  date  of  occurrence  to  the  respective  point  in  the  synodic  cycle, 
they  derived  a  frequency  distribution  that  showed  a  significant 
tendency  for  heavy  precipitation  to  occur  three  times  more 
frequently  In  the  middle  of  the  first  and  third  weeks  of  the  synodic 
month  compared  to  the  second  and  fourth  weeks.  Adderly  and  Bowen 
(1962)  found  similar  results  in  the  Southern  Hemisphere  for  Mew 
Zealand.  Vlsvanathan  (1966b)  examined  the  distribution  of  heavy 
rainfall  In  India  with  respect  to  the  synodic  period  and  found  a 
definite  lunar  component  In  the  frequency  of  occurrence  -  again  with 
two  maxima  and  minima  during  the  period.  He  also  found  a  shift  In 
phase  and  decrease  In  amplitude  of  the  precipitation  frequency  of 
occurrence  curves  with  an  Increase  In  latitude.  Bradley  (1964) 
studied  the  transition  of  269  North  Atlantic  tropical  depressions 
Into  hurricanes  and  found  the  timing  to  be  dependent  on  both  the 
half-synodic  and  the  anomalistic  (perigee  to  perigee)  periods  which 
are  the  most  Important  cycles  for  determining  the  magnitude  of  the 
Instantaneous  soli -lunar  gravitational  force.  Vlsvanathan  (1966a) 
found  that  the  frequency  of  formation  of  monsoon  depressions  In  the 
Bay  of  Bengal  and  Arabian  Sea  from  1877  to  1960  as  a  function  of  the 

lone  half  of  the  lunar  synodic  period  Is  14.7365  days  and  aliases 
at  the  frequency  .26  yr"1  In  a  spectrum  of  data  that  Is  sampled 
once  per  year. 


synodic  month  had  a  strong  14.765  day  wave  component. 

Brier  (1968)  has  also  presented  the  possibility  that 
quasi -biennial  variations  may  be  forced  by  the  soli -lunar 
gravitational  force,  for  In  an  analysis  of  the  zonal  Index  he  found 
a  significant  spectral  peak  at  3.25  years  and  was  then  able  to 
relate  the  characteristics  of  the  zonal  Index  time  series  to  those 
of  the  soil-lunar  gravitational  tidal  force. 

Thus,  there  Is  empirical  evidence  for  soil-lunar  tidal  forcing 
on  the  atmosphere.  However,  few,  if  any,  of  these  Ideas  have  found 
their  way  Into  practical  use,  mainly  because  the  arguments  are 
usually  statistical  and  linear  and  because  the  forces  Involved  are 
generally  felt  to  contribute  trivially  to  the  total  atmospheric 
response.  This  research  Is  based  on  the  hypothesis  that  soil-lunar 
tidal  forcing  does.  In  fact,  produce  an  element  of  climatic 
variability  and  will  address  a  case  where  lunar  tidal  forcing 
apparently  contributes  In  a  significant  manner  to  that  part  of  the 
variance  which  lies  In  the  Interannual  range. 

The  diurnal  and  semi-diurnal  tides  in  the  ocean  and  atmosphere 
are  well  known  manifestations  of  the  lunar  tides.  However,  as  shown 
by  Bryson  (1948)  and  discussed  by  Brier  (1968)  the  longer  period 
lunar  tides  may  be  more  Important  In  the  atmosphere  In  that  they  may 
cause  shifts  In  the  geographic  locations  of  steady  or  regular 
features,  such  as  subtropical  anticyclones.  Such  shifts  would 
result  In  large  departures  from  normal  In  the  areas  where  the 
gradients  are  the  largest.  Recently  Currie  (1981)  has  provided 


evidence  of  the  18.6  year  lunar  nodal  cycle  in  North  American 
temperature  and  drought  data. 

In  this  chapter  I  will  present  the  Information  needed  to 
understand  the  source  of  the  two  most  important  Interannual  lunar 
tidal  periodicities.  The  next  chapter  will  then  test  the  lunar 
tidal  *orc1ng  hypothesis  by  analyzing  the  response  of  station 
pressure,  precipitation,  and  the  300  mb  u-component  of  the 
geostrophlc  wind  at  the  hypothesized  forcing  frequencies  to  see  what 
relationship  exists.  If  any,  between  lunar  tidal  forcing  and  the 
atmosphere. 

2.2  Basic  tidal  analysis. 

Only  the  lunar  component  will  be  considered  here.  The  Sun's 
effect  on  tidal  forcing  will  be  discussed  when  considering  the 
variations  In  Earth-Moon  distance  In  section  2.2.2. 

Mote:  In  the  discussion  that  follows  it  will  be  necessary,  at 
times,  to  use  the  equatorial  coordinate  system.  In  this  coordinate 
system  right  ascension  and  declination  bear  the  same  relationship  to 
the  celestial  equator  and  poles  as  longitude  and  latitude, 
respectively,  bear  to  the  terrestlal  equator  and  poles.  Right 
ascension  gives  the  number  of  degrees  measured  eastward  along  the 
celestial  equator  from  the  vernal  equinox  which  Is  the  point  where 
the  celestial  equator  and  the  ecliptic  planes  Intersect. 

Declination  gives  the  number  pf  degrees  north  or  south  of  the 
celestial  equator. 


Following  the  development  given  by  Chapman  and  Llndzen  (1970), 
the  tidal  potential  due  to  the  Moon  at  a  point  P  (see  Fig  2.1)  on 
the  Earth  Is  the  difference  between  the  attraction  of  the  Moon  at 
point  P  and  the  attraction  for  the  Earth  as  a  whole.  The  local 
tidal  potential  Is: 

Up  «  GM/L  -  GM/(D2  -  2aDcosx  +  a2)1/2,  2.1 

where  G  Is  the  universal  gravitation  constant,  M  Is  the  mass' of  the 
Moon,  and  a  Is  the  radius  of  the  Earth  and  where  L,  the  distance 
from  point  P  to  the  Moon,  has  been  rewritten  using  the  Law  of 
Cosines.  The  tidal  potential  associated  with  the  Earth  Is: 

^Earth^tW02^3  cosx. 

Then  we  have: 

U tidal  *  UEarth  “  Up  *  - §M -  2.3 

D2  (D2-2aDcosx+a2)1/2 

The  first  term  of  an  expansion  In  powers  of  a/D  yields: 

Utidal  -  -3/2  GMa2(cos2x-l/3)  2.4 

IP 

which  Is  an  approximate  form  of  the  tide  generating  potential. 

The  potential  Is  then  basically  a  function  of  the  angle  between 
the  Moon  and  the  point  P.  This  angle  Is  constantly  changing  and  Is 
dependent  upon  three  basic  processes: 

1.  The  rotation  of  the  Earth  about  Its  axis.  It  takes  23  hr 
53  min  4  sec  for  the  Earth  to  rotate  once  with  respect  to  the  stars. 
This  Is  a  sidereal  day  which  Is  .9973  of  a  solar  day  (24  hrs). 

2.  The  revolution  of  the  Moon  In  Its  orbit  about  the  Earth 
from  ascending  node  to  ascending  node.  This  Is  the  nodical  month  and 


Latitudinal  geometry  used  for  the  tidal  potential 
equation  development.  The  North  Pole  Is  at  N;  the 
center  of  the  Earth  Is  at  0;  the  Moon  is  at  M;  and 
the  point  of  Interest  Is  at  P.  The  figure  Is  not 
drawn  to  scale. 


takes  27.2122  days.  The  nodes  of  the  Moon's  orbit  occur  at  the 
Intersection  of  the  Moon's  orbital  plane  and  the  ecliptic  plane. 

3.  The  retrogression  of  the  line  of  nodes  which  takes 
18.6133  years  (6798.3658  days).  Since  the  line  of  nodes  Is 
retrogressl ng  at  the  same  time  the  Moon  Is  revolving  about  the  Earth, 
when  the  Moon  Is  back  to  Its  ascending  node  again  It  Is  still  not  back 
to  the  same  place  In  Its  orbit  with  respect  to  the  stars.  This 
period,  called  the  sidereal  month,  takes  slightly  longer  than  the 
nodical  month  and  Is  27.3216  days  long. 

Let  a  and  e  be  the  colatitude  of  the  sublunar  point  and  the 
location  of  Interest,  respectively,  as  shown  In  Figure  2.1.  Then  If 
we  let  i  be  the  east  longitude  of  point  P  and  o  be  the  hour  angle  of 
the  Moon  to  the  west  of  the  Greenwich  meridian  (see  Fig.  2.2),  the 
cos*  can  be  written  as: 

cosx«cosAcose+s1nAs1necos(a+0).  2.5 

The  angles  e  and  i  are  related  to  point  P's  latitude  and  longitude, 
respectively.  The  angle,  a,  can  be  represented  as  follows: 

a  «  w/2-(23.45*+5.15*cosNo)cosSm,  2.6 

where  23.45*  Is  the  latitude  of  Intersection  of  the  ecliptic  plane 
with  the  earth, 

5.15*  Is  the  angle  at  which  the  Moon's  orbital  plane  Is  Inclined 
to  the  ecliptic  plane. 

No  Is  2rf,^t  where  fu_  Is  the  nodal  frequency  of  .0537 


Moon 


ine  of  Nodes 


Figure  2.2.  Longitudinal  geometry  used  for  the  tidal  potential 
equation  development.  The  outer  circle  represents 
the  Moon's  orbit  and  the  inner  circle  represents  the 
Earth  with  center  at  0.  P  is  the  point  of  Interest. 
The  figure  is  not  drawn  to  scale. 


Sm  Is  2*fSmt  where  Is  the  sidereal  month  frequency 
13.3683  year'll. e.,  365.2422/27.3216). 

Next  let's  consider  the  hour  angle,  a.  If  the  Moon  were 
stationary,  then  the  only  thing  affecting  a  would  be  the  rotation  of 
the  earth  and  a  would  go  through  360*  In  one  sidereal  day.  In  this 
case  a«2»t/.9973  with  t  measured  In  solar  days.  However,  the  Moon 
revolves  around  the  Earth  as  shown  In  Figure  2.2.  It  takes  the  Moon 
27.2122  days  to  travel  west  to  east  around  the  Earth  from  ascending 
node  to  ascending  node,  so  If  one  were  to  measure  a  at  the  end  of  each 
Earth  rotation  It  would  be  smaller  In  terms  of  degrees  west  of  the 
Greenwich  meridian.  At  this  point, 

a»(2»/. 9973  -  2*/27.2122)t.  2.7.1 

There  Is  still  one  other  factor  to  consider  and  that  Is  the 
retrogression  of  the  line  of  nodes  from  east  to  west  as  depicted  in 
Figure  2.2.  This  retrogression  Increases  a  so  we  have, 

o-(2ir/.9973  -  2t/27.2122  ♦  2w/6798.4)t.  2.7.2 

Alpha  can  also  be  represented  as: 

a*>(2*/.9973  -  2^/27.255241  -  2w/aps1des)t,  2.8 

where  Instead  of  the  nodical  month  and  the  nodal  period  the 
anomalistic  month  and  the  period  of  the  revolution  of  the  line  of 
apsides  are  used.  The  line  of  apsides  joins  the  points  of  perigee  and 
apogee  and  makes  one  revolution  In  8.849  years. 

Now  In  order  to  be  able  to  calculate  the  potential  at  a  given 
location  the  expression  (cos  *  -  1/3)  must  be  evaluated.  It  Is  the 
sum  of  the  following  three  terms: 
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Term  1:  (3/2) (cos2a  -  l/3)(cos2e  -  1/3)  2.9 

Term  2:  (l/2)s1n2As1n2ecos(a+d)  2.10 

Term  3:  (l/2)s1n2AS1n2ecos2(a+d)  2.11 

2.2.1  Calculation  of  tidal  potential. 

Terms  1,  2,  and  3  have  been  evaluated  to  get  a  physical  Idea  of 
what  Influences  the  value  of  the  lunar  potential.  Certain  simplifying 
assumptions  have  been  made  In  the  calculations: 

1.  The  distance  between  Moon  and  Earth  Is  constant  (l.e. ,  D 
In  Eqn.  2.4  Is  constant).  See  section  2.2.2  for  a  discussion  of  the 
varying  Earth-Moon  distance. 

2.  In  Eqn  2.6  the  value  of  (23.45°  +  5.15°cosNo)  was  assumed 
constant  for  a  given  month  and  year.  The  maximum  declination  of  the 
Moon  goes  through  an  18.6  year  cycle,  thus  the  value  Is  fairly 
constant  for  a  30  day  period.  This  value  was  taken  to  be  the  maximum 
northward  declination  for  the  month  and  year  with  the  phasing  on  the 
cosSm  term  of  Eqn.  2.6  being  the  time  of  this  maximum  northward 
declination.  Both  the  declination  and  the  time  were  found  In  the 
American  Ephemerfs  and  Nautical  Almanac  (1855-1982). 

In  order  to  calculate  the  hour  angle,  a,  the  following  equations 
taken  from  the  American  Ephemerls  and  Nautical  Almanac  (1980)  were 
used: 


P'«  281 *.22083  ♦  0.0000470684d  +  0.000453T2  +  0.000003T3  2.15 
N  »  259*. 183275  -  0.0529539222d  +  0.002078T2  +  0.000002T3  2.16 
e  -  0.054900489  2.17 
1  «  5*. 145396  2.18 
where 


L  Is  the  mean  longitude  of  the  Moon,  measured  In  the  ecliptic  from 

the  vernal  equinox  to  the  ascending  node  and  then  along  the 

orbit  to  the  Moon. 

P  Is  the  longitude  of  perigee  of  the  Moon's  orbit. 

N  Is  the  longitude  of  the  Moon's  ascending  node. 

1  Is  the  Inclination  of  the  Moon's  orbit  to  the  ecliptic. 

e  Is  the  average  eccentricity  of  the  Moon's  orbit. 

L' ,  P*  are  similar  to  L  and  P  except  they  pertain  to  the  Sun. 

In  these  equations  T  Is  measured  in  Julian  centuries  from  1900  January 
0.5  Ephemerls  Time  and  d  Is  the  number  of  ephemerls  d^ys  from  this 
date.  (An  ephemerls  day  Is  86400  ephemerls  seconds  where  an  ephemerls 
second  Is  1/31556925.9747  of  the  tropical  year.) 

These  formulas  must  then  be  corrected  for  certain  well  known 
perturbations  of  the  Moon's  orbit.  Voigt  (1974,  p.  59)  discusses 
evectlon,  variation  and  the  annual  equation.  Goldstlne  (1973)  made 
use  of  the  following  equations  to  take  these  perturbations  into 
account. 

Ax  .  4467*s1n(L  -  2L '  +  P)  2.19 

Aj  .  2145"s1n2(L  -  L')  2.20 

A3  .  658"s1n(L'  -  P'  ♦  180*) 


2.21 


2.22 


A4  .  198"s1n(L  -  3L*  ♦  P  +  P') 

Ag  .  155"s1n(2L  -  3L'  +  P')  2.23 

L  can  then  be  replaced  by  the  corrected  mean  longitude: 

LM  «  L  ♦  Aj  ♦  A2  +  A3  ♦  A4  ♦  Ag.  2.24 

In  order  to  find  the  true  longitude,  y»  of  the  Moon  In  Its  orbit 
the  equation  of  center,  as  discussed  by  Voigt  (op.clt.,  p .59)  and 
Goldstlne  (op.  clt. )  Is  used.  When  retaining  terms  through  e4,  the 
equation  becomes: 

y  »  L  ♦  (180/ ■»)  [(2e  -  e3/4)s1n(L-P)  + 

(5e2/4-lle4/24)s1n2(L-P)  +  (13e2/12)s1n3(L-P) 
(103e4/96)s1n4(L-P)].  2.25 

Thus,  we  now  have  the  true  longitude  of  the  Moon's  orbit.  This 
angle  Is  measured  from  the  vernal  equinox  to  the  ascending  node  In  the 
plane  of  the  ecliptic  and  then  along  the  orbit  to  the  Moon's  true 
position.  The  Moon's  position  must  then  be  projected  onto  the  plane 
of  the  ecliptic  using  the  approximate  formula  known  as  the  reduction 
to  the  ecliptic  found  In  Brower  and  Clemence  (1961,  p.  47): 

SL  .  L  ♦  (180/t)  [p2s1n2(Y  -  N)  -  (p4/2)s1n4(Y  -  N) 

♦  (p6/3)s1n6(T  -  N)]  2.26 

where  p«tan(1/2) 

All  of  the  above  equations  were  used  to  find  the  change  In  the 
angle,  a,  during  the  course  of  the  month.  This  change  will  be 
denoted  as  ao^.  Thus,  by  knowing  the  hour  angle  at  the  first 
Instant  of  each  month  (oq),  the  sidereal  day,  and  Aat,  the  hour 
angle  at  any  time  during  the  month  can  be  calculated  by: 


o  «  a0  +  (360*/.9973)t  -  Ac»t.  2.27 

The  declination  of  the  Moon,  which  Is  the  number  of  degrees  north 
or  south  of  the  celestial  equator,  was  calculated  using  the  time  of 
maximum  northward  declination  during  the  month  and  the  value  of  the 
declination  at  that  time.  Thus,  to  find  the  colatitude  of  the 
declination,  a,  at  any  given  time  during  the  month  the  following 
equation  was  used: 

a  »  90.  -  «  «  90.  -  «mavcos[(2ir/Sm)(t-d  v)]  2.28 

niaA  IflaX 

where 

*max  «*x1mum  declination  In  days. 

6  is  the  declination  of  the  Moon, 

amax  the  max,|limim  northward  declination  at  time  dnax* 

Sm  Is  the  sidereal  month,  27.3216  days,  and 
t  Is  the  time  of  the  month  In  days. 

With  the  values  for  a  and  a  we  are  then  able  to  calculate  the  average 
potential  for  a  given  month. 

The  values  for  <xq,  and  amax  were  obtained  from  the 

American  Ephemerls  and  Nautical  Almanac  (1855-1982). 

It  Is  Important  to  emphasize  that  we  are  calculating  the  monthly 
average  potential  because  the  climatic  data  we  will  be  analyzing  are 
all  monthly  averages.  This  makes  the  external  forcing  data  set 
comparable  to  the  observed  data. 

2.2.2  Effect  of  varying  Earth-Moon  distance  on  the  forcing. 

In  the  preceding  analysis  the  distance  between  Moon  and  Earth  was 


assumed  to  be  constant,  producing  a  constant  amplitude.  In  fact, 
because  of  the  eccentricity  of  the  Moon's  elliptical  orbit  (e=.0549) 
around  the  Earth  (with  the  Earth  at  one  focus),  the  Moon  has  a  closest 
approach  to  the  Earth  during  the  month  called  perigee  and  about  two 
weeks  later  It  reaches  Its  furthest  distance  called  apogee.  The  time 
from  perigee  to  perigee  Is  known  as  the  anomalistic  month  of  27.55455 
days.  In  addition,  the  orbital  axis  joining  perigee  and  apogee,  the 
line  of  apsides,  rotates  In  a  counterclockwise  direction  completing 
one  revolution  In  8.849  years.  Because  of  the  rotation  of  the  line  of 
apsides,  the  point  In  the  lunar  orbit  where  perigee  occurs  is 
constantly  changing. 

Looking  at  the  tidal  potential  (Eqn.  2.4  and  2.9  -  2.11),  one  can 
see  that  a  variation  In  distance  will  cause  the  amplitude  of  the 
periodic  cosine  terms  to  change  significantly  In  instantaneous  value 
since  the  amplitude  Is  a  function  of  the  cube  of  the  distance. 

When  the  Moon  Is  at  perigee.  It  Is  closest  to  the  Earth  in  its 
orbit  and,  thus,  the  tide  raising  forces  are  the  strongest.  This 
effect  can  be  enhanced  further  If  perigee  occurs  at  the  time  of 
Moon-Earth-Sun  alignment  known  as  sytygy  (l.e.,  full  or  new  moon)  for 
then  the  gravitational  attractions  of  the  loon  and  Sun  act  to 
reinforce  each  other.  When  this  close  alignment  of  perigee  and  sy?ygy 
occurs  (l.e.,  the  line  of  apsides  and  line  of  sy?ygies  are  nearly 
coincident)  a  perturbation  called  lunar  evection  Increases  the 
eccentricity  of  the  lunar  orbit  which  In  turn  decreases  the  Earth-Moon 
distance  at  perigee  with  a  subsequent  increase  In  the  gravitational 


attraction.  As  a  result,  the  tidal  potential  at  the  sublunar  point 
will  be  even  more  negative  than  usual.  However,  since  the 
eccentricity  of  the  Moon's  orbit  Is  Increased,  It  not  only  makes 
perigee  closer  but  It  also  makes  the  point  of  apogee  farther  from  the 
Earth  producing  an  opposite  effect  on  the  monthly  average  lunar 
potential. 

Given  a  perlgee-syzygy  event  there  are  basically  three  factors 
which  can  further  enhance  this  effect:  1.)  the  phase  of  the  Moon 
(l.e.,  full  or  new),  2.)  the  dec! i natl onal  planes  of  the  Moon  and  Sun, 
and  3.)  whether  the  Earth  Is  at  perihelion.  The  phase  of  the  Moon  Is 
Important  because  the  gravitational  attraction  of  the  Sun  on  the  Moon 
Is  not  the  same  at  full  and  new  moon.  At  full  moon  the  gravitational 
attraction  of  both  the  Earth  and  Sun  on  the  Moon  are  in  the  same 
direction  and  causes  the  Moon  to  move  to  a  position  closer  to  the 
Earth  than  If  just  the  Earth  were  present.  At  new  moon  the 
gravitational  attraction  of  the  Earth  and  the  Sun  on  the  Moon  are 
directed  In  opposite  directions  plus  the  Moon  Is  now  closer  to  the  Sun 
than  at  full  moon.  As  a  result  the  lunar  evectlon  perturbation  Is  not 
as  great  and  the  distance  at  perigee  Is  not  decreased  as  much  as  It  Is 
with  full  moon  at  perigee. 

The  only  requirement  for  syzygy  to  occur  Is  that  the  Moon  and  Sun 
be  In  the  same  celestial  longitude.  However,  If  they  also  are  In  the 
same  dec! 1 natl onal  plane  the  combined  gravitational  attraction  of  the 
Earth  and  Sun  on  the  Moon  Is  noticeably  Increased  at  full  moon, 
Increasing  the  eccentricity  of  the  orbit  and  bringing  the  Moon  even 


closer  to  the  Earth. 

The  Importance  of  perihelion  Is  that  the  above  Influences  are  all 
enhanced  when  the  Sun  Is  closest  to  the  Earth  and  thus  to  the  Moon. 

In  fact.  Wood  (1976)  has  tabulated  all  of  the  extreme  perfgee-sysygy 
events  over  the  400  year  period  from  1600  through  1999  and  found  that 
they  all  occurred  near  perihelion  (Oct  31  through  Mar  8),  at  full  moon 
when  the  Sun  and  Moon  were  nearly  In  the  same  decllnatlonal  plane. 
Voigt  (1974)  has  stated  that  the  minimum  Earth-Moon  distance  occurs 
under  the  following  conditions: 

a.  Full  Moon 

b.  Perigee 

c.  Moon  between  the  nodes 

d.  Earth  at  perihelion 

At  perihelion  the  Sun  Is  at  about  281*  In  right  ascension  so  the  full 
moon  would  be  at  101*  which  puts  It  near  Its  maximum  positive 
declination  and  since  the  Sun  Is  at  about  Its  maximum  negative 
declination,  having  just  passed  through  the  winter  solstice,  the  two 
bodies  are  nearly  coplanar. 

Having  examined  the  main  factors  Influencing  the  Earth-Moon 
distance.  It  remains  to  be  seen  what  effect  this  has  on  the  forcing. 

If  we  now  define  this  distance,  D,  as  the  average  distance  plus  some 
sinusoidal  variation  (e.g.,  the  anomalistic  month)  then: 

D3  .  (5  ♦  D'cosA)3 

.  D3  ♦  3D2D'cosA  ♦  3D,25cos2A  ♦  D,3cos3A 
.  D3[l+3(O'/0)cosA  ♦  3(D,2/02)cos2A 
♦  (D,3/03)cos3A] 


2.29 


2.30 


Letting  h«D'/D  and  x^cosA  we  have 

D-3«D~3[1  ♦  3hx  ♦  3h2x2  +  h3x3]-1. 

Using  a  Taylor  Series  expansion  of  D  about  0  yields: 

D_3»D"3[1  -  3hx  ♦  6h2x2  +  ....]+1  2.31 

By  retaining  only  the  second  order  terms  equation  2.4  becomes: 

Ut1da1““(3/2> (G"»a2/D3)C1  -  (3D'/6)cosA  ♦  2.32 

(6D,2/D2)cos2A][cos2x  -1/3]. 

Ue  will  see  that  the  nodal  frequency  dominates  the  forcing. 

Looking  then  at  the  varying  part  of  equation  2.32  we  have: 

U* tidai  -  [1-(3D'/6)cosA  +  (6D,2/D2)cos2A]cosNo.  2.33 

The  maximum  Instantaneous  distance  variation  possible  Is  27,990  km. 
Thus  D'/D  Is  .0728  and  this  value  squared  Is  .0053  making  the  main 
contribution  to  the  unmodulated  nodal  cycle.  When  dealing 

with  monthly  averages,  the  D'  value  Is  much  less  making  the  effect  of 
distance  variation  seemingly  minimal  on  the  monthly  average  of  the 
forcing.  Thus,  even  though  D  was  a  constant  In  the  calculations,  this 
did  not  have  any  great  Impact  on  the  monthly  average  potential  values. 

2.2.3  Magnitude  of  lunar  tidal  potential. 

At  a  given  location  on  the  Earth,  Term  1  (Eqn.  2.9)  Is  only 
dependent  on  the  declination  of  the  Moon,  while  Terms  2  and  3  are  also 
dependent  upon  the  longitude  of  the  sublunar  point.  The  Moon  passes 
over  the  same  longitude  once  every  24  hr  50.47  min.  Thus,  over  that 
period  the  part  of  the  potential  determined  by  Terms  2  and  3  nearly 
averages  out  to  zero.  It  does  not  average  out  exactly  to  zero. 


however,  since  the  latitude  of  the  sublunar  point  changes  slightly 
over  that  time.  This  can  be  seen  by  looking  at  Table  2.1  which  shows 
the  Instantaneous  values  of  the  three  terms  at  various  configurations 
of  the  Moon's  declination  and  hour  angle.  The  global  pattern 
corresponding  closely  to  Table  2.1.1  Is  seen  In  Figure  2.3  which  Is 
the  Instantaneous  potential  for  August  0.0  1948.  The  sublunar  point 
Is  at  24.85*M,  239.67*W.  The  slight  difference  between  these  numbers 
and  those  In  Table  2.1.1  exists  because  the  table  assumes  a  constant 
Earth-Moon  distance  whereas  the  figure  uses  Doodson  (1921)  harmonics 
which  have  the  effect  of  the  variable  Earth-Moon  distance  Included. 

Instantaneously,  Terms  2  and  3  can  be  an  order  of  magnitude  larger 
than  Term  1  as  seen  In  Table  2.1.  However,  the  magnitude  of  Terms  2 
and  3  are  two  orders  of  magnitude  less  than  Term  1  when  a  monthly 
average  Is  considered.  For  example,  the  values  for  June  1904  were: 
Term  1,  -0.34174  m2s“2;  Term  2,  -0.00165  m2s~2;  and  Term  3, 

0.00049  m2s*2. 

One  Interesting  point  to  note  about  Term  1  Is  that  at  35.26*N  and 
at  35.26*$  the  term  goes  to  zero.  Poleward  of  this  latitude  the 
factor,  (cos  e  -  1/3)  of  Eqn.  2.9,  Is  positive  and  equatorward  It  Is 
negative.  Terms  2  and  3  are  zero  at  the  poles  while  Term  2  Is  also 
zero  at  the  equator. 

2.2.4  Affect  of  geographic  location  on  each  term. 

2. 2. 4.1  Term  1. 

A  look  at  Equation  2.9  for  Term  1  shows  that  it  has  no  longitude 


Figure  2.3.  Instantaneous  lunar  potential  for  the  first  Instant  of  August  1948 
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dependence  and  that  It  varies  latltudlnally  by  the  factor, 

(cos  e-1/3),  where  e  Is  the  colatitude  of  the  location.  Thus  all 
points  between  35.26*S  and  35.26*N  are  In  phase  with  each  other  and 
180*  out  of  phase  with  all  points  poleward  of  these  latitudes. 

2. 2. 4. 2  Term  2  and  Term  3. 

Both  of  these  terms  have  a  latitudinal  and  longitudinal 
dependence.  Looking  at  the  factor  cos{a+d)  and  considering  the  hour 
angle  a  as  the  sum  of  the  Initial  hour  angle  at  the  beginning  of  the 
month,  oq,  plus  the  change,  c^,  up  to  time  t  we  have 
cos(at+a0+d).  Thus,  the  phase  for  this  term  is  OQ+d.  In  other 
words,  the  timing  of  maximum  and  minimum  values  at  a  location  Is 
dependent  on  both  the  Initial  hour  angle,  aQ,  and  the  longitude,  d, 
of  the  point. 

2.2.5  Important  frequencies  associated  with  each  term. 

The  calculations  described  above  were  performed  for  a  location  at 
25*N,  75*E  for  the  years  1899  through  1960  for  June.  Five  monthly 
average  time  series  were  generated: 

Total  lunar  potential  U  (»Uj  ♦  U2  +  U3) 

Term  1  potential*  U* 

Term  2  potential,  U2 
Term  3  potential,  Uj 
Term  2  ♦  Term  3  potential,  U23* 

To  see  which  frequencies  were  predominant  a  perlodogram  of  each  time 
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series  was  constructed.  The  results  are  summarized  In  Table  2.2. 

Note:  The  perlodogram  was  calculated  using  tenth  band  Intervals, 
therefore,  with  a  62  year  time  period  a  resolution  of  .0016  year'1 
was  attainable.  Thus,  for  the  frequency  .0532  year-1,  the  next 
lower  frequency  would  be  .0516  year"1  and  the  next  higher  frequency 
would  be  .0548  year  -1. 

This  table  shows  that  the  two  frequencies  .0532  year"1  and  .2645 
year'1  explained  almost  100  percent  of  the  variance  for  U.  For 
the  two  frequencies  were  .0532  year"1  and  .2629  year"1.  Most  of 
the  variance  for  Ug  and  Uj  Is  explained  by  the  frequency  between 
.2629  year"1  and  .2645  year'1.  Thus,  all  three  terms  exhibit  a 
frequency  around  .263  year"1;  however,  looking  at  the  amplitudes  of 
this  frequency  for  the  three  terms,  we  see  that  Uj's  amplitude  Is 
over  three  times  that  of  UZ3. 

The  origin  of  the  two  frequencies  .0537  year"1  and  .263  year'1 
will  be  examined  In  the  next  two  sections. 
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2.3  Lunar  tidal  potential  at  the  nodal  frequency  (.053725  year"1) 

As  Table  2.2  shows,  the  nodal  cycle  Is  the  predominant  frequency 
In  the  monthly  mean  lunar  potential  time  series  and  explains  over  95 
percent  of  the  Interannual  variance  for  the  month  of  June. 

The  nodal  cycle  of  18.6  year  period  Is  the  time  It  takes  the  pole 
of  the  Moon's  orbit  to  rotate  about  that  of  the  ecliptic.  The  Moon's 
orbital  plane  Is  Inclined  to  the  ecliptic  at  an  average  angle  of 
5a9'.  Figure  2.4  taken  from  Wood  (1976)  shows  how  the  nodal  cycle 
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TABLE  2.2 


Important  Frequencies  for  June  Lunar  Potential  at  25#M,  75#E 


Term 

F^a 

Percent  Var/ 
Amplitude 

U 

.0532 

96.9/4.0E-02 

«1 

.0532 

95. 5/4. OE -02 

U2 

— 

U3 

— 

— 

u23 

-- 

-- 

Percent  Var/ 


Phase 

Freq 

Ampl 1 tude 

96.9 

.2645 

2. 7/6. 6E -03 

96.3 

.2629 

4.0/8. IE-03 

.2629 

89.6/2. 7E-03 

— 

.2645 

76. 9/8. OE -04 

— 

.2645 

85.1/2.5E-03 

PI 

1< 

-II 

l 
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varies  the  latitudinal  range  over  which  the  Moon  traverses  during  one 
sidereal  month.  His  diagram  assumes  an  Inclination  of  5°  for 
simplicity. 

When  the  ascending  node  Is  at  the  vernal  equinox  (l.e.,  the 
subsolar  point  on  the  celestial  sphere  at  the  time  of  the  vernal 
equinox),  the  Moon  achieves  Its  maximum  declination  of  28.6°  north  and 
south;  whereas  9.3  years  later  when  the  ascending  node  Is  at  180°  from 
the  vernal  equinox,  the  maximum  lunar  declination  Is  18.3°  north  and 
south. 

As  Table  2.1  and  Figure  2.3  show,  the  Instantaneous  potential  at 
the  sublunar  point  Is  negative.  If  during  the  course  of  the  sidereal 
month,  the  Moon  Is  constrained  to  the  ^  18.3*  latitudinal  range  then 
the  average  potential  for  the  month  will  be  lower  In  this  region  than 
If  It  had  a  larger  range,  for  the  large  area  of  negative  potential 
seen  In  Figure  2.3  Is  also  constrained.  Table  2.3  shows  this  effect. 
In  the  table  are  the  calculated  values  of  the  average  June  lunar 
potential  for  the  years  1903-1922  at  25*N,  75*E  along  with  the  maximum 
declination  for  that  particular  June.  As  one  can  see  the  further 
north  the  sublunar  point  can  be  the  less  negative  the  average 
potential.  The  lack  of  an  exact  correspondence  between  the  average 
potential  and  the  maximum  declination  Is  because  of  an  effect  related 
to  the  date  during  the  month  on  which  the  Moon  reaches  Its  maximum 
declination.  This  will  be  discussed  In  section  2.4. 

Since  the  nodal  cycle  arises  from  Term  1,  this  cycle  has  no 
longitude  dependence.  Thus,  In  years  when  the  ascending  node  Is  near 


EFFECT  OF  MOON’S  18.6-YEAR  NODICAL  CYCLE 
UPON  THE  MAXIMUM  LUNAR  DECLINATION 


Effect  of  the  nodal  cycle  on  the  latitudinal  range  of  the  sublunar  point 
(from  Wood,  1976). 


TABLE  2.3 


Relationship  of  Average  Potential  to  «max  and 

JUNE  6mav  TI$ia&F 

AVERAGE  POTENTIAL  MAX  DECLINATION  MAX  DECLINATION 

1903  -.3478  m2s‘2  18.58#N  day  23.38 


the  vernal  equinox  the  latitudinal  distribution  of  the  lunar  potential 
anomaly  due  to  the  nodal  cycle  Is  as  seen  In  Fig.  2.5.1  and  just  the 
opposite  when  the  ascending  node  Is  at  the  autumnal  equinox  as  In 
Figure  2.5.2. 

2.4  Lunar  tidal  potential  at  frequency  of  .263  year"1. 

Once  the  nodal  cycle  has  been  considered,  the  only  other  frequency 
that  explains  more  than  one  percent  of  the  variance  Isa  frequency 
around  .263  year"1-  There  are  two  parameters  of  the  Moon's  orbit 
for  a  given  month  that  determine  the  lunar  tidal  potential  anomaly  at 
this  frequency.  They  are  the  time  of  maximum  declination  and  the  hour 
angle  at  the  first  Instant  of  the  month. 

2.4.1  Time  of  maximum  declination,  dmax. 

2.4. 1.1  Influence  on  Term  1. 

Table  2.3  shows  the  day  Umax)  when  the  Moon  Is  at  Its  maximum 
declination  for  the  month  of  June  In  the  years  Indicated.  This  timing 
Is  the  single  most  Important  factor  In  determining  the  magnitude  of 
the  lunar  potential  monthly  average  for  Term  1  If  the  nodal  cycle 
Influence  Is  held  constant.  As  already  discussed,  Term  1  Is  the 
largest  of  the  three  terms  by  a  factor  of  100  when  monthly  averages 
are  calculated  and  the  amplitude  of  this  frequency  Is  over  three  times 
that  of  the  other  two  terms  combined.  Therefore,  by  considering  just 
this  term  we  can  see  the  effect  that  the  time  of  maximum  declination 
has  on  the  tidal  potential. 


POTENTIAL  ANOMALY  n  POTENTIAL  ANOMALY 


FIGURE  2-S.l  -  NOOE  AT  0  DEGREES  RIGHT  ASCENSION 


IGURE  2-5.2  -  NOOE  AT  180  DEGREES  RIGHT  ASCENSION 


Figure  2.5.  Latitudinal  distribution  of  the  lunar  potential  anomaly 
at  the  extremes  of  the  nodal  cycle  (units  are 
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Term  1  Is  dependent  on  the  declination  of  the  Moon  (a  In  Fig  2.1 
and  Eqn.  2.9)  only  and  not  on  the  hour  angle.  Since  a  Is  a  function 
of  the  cosine  of  the  sidereal  month  (see  Eqn.  2.28),  the  average  value 
of  a,  for  a  given  <5max,  after  27.3216  days  Is  the  same  no  matter 
what  value  dmax  has.  Thus,  It  Is  the  remaining  days  of  the  month 
that  determine  If  the  potential  will  be  above  or  -low  the  first 
27.3216  day  average.  If  the  Moon  begins  the  month  lear  maximum  north 
declination.  It  will  end  the  month  In  that  position  also.  This 
configuration  produces  the  maximum  potential  (smallest  negative). 
Likewise  If  the  Moon  starts  the  month  at  maximum  south  declination, 
the  same  result  Is  achieved  since  the  tidal  potential  at  the  sublunar 
point  Is  nearly  matched  on  the  opposite  side  of  the  earth  (e.g.  see 
Figure  2.3).  The  minimum  monthly  average  potential  occurs  when  the 
average  latitude  of  the  Moon  for  the  remaining  days  of  the  month  Is 
over  the  equator.  The  reason  for  this  can  be  seen  qualitatively  from 
Figure  2.3.  In  this  figure  the  sublunar  latitude  Is  about  25*N.  To 
calculate  an  approximate  lunar  potential  average  for  a  location  at 
latitude  25*M  for  the  day  depicted  In  this  figure,  all  one  need  do  Is 
average  the  Instantaneous  values  across  all  longitudes.  This  average 
Is  a  smaller  negative  number  than  one  would  obtain  If  the  Moon  were 
over  the  equator.  In  this  case  both  large  areas  of  negative  potential 
In  Figure  2.3  would  be  displaced  toward  the  equator  producing  a  more 
negative  dally  average  at  25*N. 

Since  the  number  of  days  Involved  In  the  residual  beyond  the 
sidereal  month  varies  depending  on  the  length  of  the  month,  the  time 
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of  maximum  declination  that  produces  the  maximum/minimum  potential 
varies  slightly  from  month  to  month.  For  a  month  with  30  days  the 
following  values  for  the  time  of  maximum  declination,  dmax,  pertain 
for  a  location  between  35.26*N  and  35.26*S: 

maximum  average  monthly  potential  with  dmax  either  at  1.34  or 
15.00  days,  and 

minimum  average  monthly  potential  with  d(nax  either  at  8.17  or 
21.83  days.  Thus,  the  largest  of  the  three  terms  contributing 
to  the  total  potential  Is  very  dependent  on  the  timing  of  maximum 
declination  during  a  given  month. 

2.4. 1.2  Origin  of  the  .263  year"1  frequency. 

A  perlodogram  of  the  time  series  (1899-1960)  of  for  the 
month  of  June  yields  the  two  largest  peaks  at  .2629  year"1  (15.2 
percent  of  the  variance)  and  at  .3677  year"1  (61.4  percent). 

Nonlinear  regression  more  accurately  defines  these  peaks  at  .26308 
year"1  and  .36834  year’1. 

The  equation  relating  the  declination  and  right  ascension  of  the 
Moon  Is: 

sins  »  0.406s1na  ♦  0.008s1n3a  +  0.090s1n(a-N)  +  0.006s1n(3a-N), 

2.34 

where  «  Is  the  declination, 

a  Is  the  right  ascension, 

N  Is  the  longitude  of  the  ascending  node. 

This  shows  that,  by  considering  the  first  and  by  far  the  largest  term 
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on  the  right  hand  side,  the  Moon  attains  Its  maximum  northward 
declination  at  a  right  ascension  of  90*  and  Its  maximum  southward 
declination  at  a  right  ascension  of  270*.  Thus,  the  point  of  Its 
maximum  northward  declination  Is  always  at  the  same  point  In  the  orbit 
with  respect  to  the  stars.  This  period,  from  90*  back  to  90*  again. 

Is  the  sidereal  month  whose  value  Is  determined  by  considering  the 
nodical  month  (Mm)  and  nodal  cycle  (No)  as  follows: 

Sidereal  lunar  month  ■  (1/Nm  -  1/No)'1  2.35 

.  (1/27.21222  -  1/6798. 36575) -1  «  27.321582  days 
In  one  tropical  year  (365.2422  days)  the  Moon  Is  thus  at  Its  maximum 
northward  declination  13.3683  times.  By  sampling  this  cycle  only  once 
every  year  we  see  the  alias  at  .3683  year-1  as  was  found  In  the 
perlodogran  and  nonlinear  regression  of  the  dmax  time  series. 

In  actuality  the  Important  frequency  Is  the  first  harmonic  of 
.3683  year'1  since  the  spectral  analysis  does  not  take  Into  account 
that  the  same  average  potential  results  when  the  Moon  Is  at  Its 
maximum  southward  declination  as  discussed  In  the  previous  section. 
There  we  saw  that  a  maximum  northward  declination  time  of  about  1.34 
days  yields  the  same  average  monthly  potential  as  a  maximum 
declination  time  of  15.0  days.  We  are  now  looking  at  a  period  of 
13.6608  (27.321582/2)  days  which  goes  through  26.736534  cycles  during 
a  tropical  year.  This  frequency  aliases  at  .2635  year-1  when 
sampling  Is  done  once  a  year.  This  then  Is  the  frequency  to  expect  to 
find  In  a  time  series  of  potential  and  can  be  found  by  calculating  the 
first  harmonic  of  the  peak  near  .368  year'1. 


In  the  case  of  June 


(1899-1960)  the  peak  Is  at  .36834  year"1  giving  a  first  harmonic  at 
.7367  year"1  which  aliases  at  .2633  year"1. 

For  August  the  two  largest  peaks  In  the  perlodogram  of  are 
at  .2662  year"1  (14.1  percent)  and  at  .3677  year"1  (57.9 
percent).  Refinement  by  nonlinear  regression  yields  .2655  year"1 
and  .3674  year  -1.  The  first  harmonic  of  the  largest  peak  aliases 
at  .2652  year"1.  Thus,  there  Isa  difference  In  the  frequency  to 
expect  for  June  and  August. 

The  relative  latitudinal  distribution  of  the  anomalies  at  this 
frequency  will  be  the  same  as  In  Figure  2.5  because  both  frequencies 
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cone  from  Term  1  which  has  a  (cos  e  -  1/3)  latitude  dependence.  The 
only  difference  Is  that  the  .263  year"1  anomalies  have  a  smaller 
magnitude. 

2.4.2  Hour  angle  of  the  Moon  at  the  beginning  of  the  month,  aQ. 

Figure  2.6  shows  a  scatter  diagram  of  hour  angle,  aQ,  versus 

for  the  months  of  June  and  August.  The  data  were  taken  from 
max 

the  American  Ephemerls  and  Nautical  Almanac.  In  both  months  It  Is 
clear  that  associated  with  each  dmax  there  Is  an  oQ;  however,  the 
associated  oq  Is  different  from  one  month  to  another.  Between  June 
and  August  the  difference  Is  60  to  70  degrees.  This  difference  pleys 
an  important  role  In  Terms  2  and  3  which  are  dependent  on  the  hour 
angle  of  the  Moon. 

As  we  have  seen  In  Table  2.2  the  frequency  around  .26  year"1 
appears  In  all  three  terms.  We  have  now  also  seen  that  0^,  which 


Figure  2.6.  Scatter  diagram  of  the  hour  angle  of  the  Moon  at  the 
beginning  of  the  month  versus  the  day  of  maximum 
declination. 
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determines  the  .26  year"*  phase  In  Term  1,  Is  associated  with  oq 
which  In  turn  affects  when  Terms  2  and  3  win  be  a  maximum  at  a  given 
longitude.  Thus,  If  «.q  produces  a  maximum  value  at  the  frequency 
.26  year-*  In  Terms  2  and/or  3  at  the  same  time  that  |L.V  produces 
a  maximum  at  that  frequency,  then  the  total  potential  will  have  a 
stronger  si gnal  at  .26  year"* . 

2.5  Component  Force  (Fx  and  Fy)  derived  from  lunar  potential. 

Taking  the  negative  gradient  of  the  lunar  potential  field  results 
In  the  associated  horizontal  tidal  force: 

F  ■  -dU/dS  where  S  is  along  the  gradient. 

The  tidal  force  anywhere  on  the  Earth  can  be  resolved  Into  a  component 
vertical  to  the  Earth's  surface  and  another  horizontal  to  It.  The 
horizontal,  or  tractive,  force  Is  the  one  that  produces  the  tides.  We 
shall  briefly  consider  both  the  X  and  Y  components  of  the  horizontal 
force. 

2.5.1  F  component. 

Term  1  Is  the  major  contributor  to  this  component  of  the  force. 

As  noted  earlier  there  Is  a  nodal  point  In  the  potential  at  35.26* 
north  and  south.  This  establishes  the  variability  of  the  potential  as 
a  standing  wave.  As  seen  before  this  variability  Is  dominated  by  two 
frequencies,  .053725  year"*  and  .263  year"*.  Figure  2.5  showed 
the  latitudinal  distribution  of  the  potential  at  the  extremes  of  the 
nodal  cycle.  Now,  since  the  force  Is  the  negative  gradient  of  the 
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potential  we  see  the  configuration  of  Figure  2.5.1  gives  a  positive 
Y-component  (l.e.,  a  northward  directed  force)  and  the  opposite 
configuration  would  give  a  negative  Y-component  (l.e,  a  southward 
directed  force). 

The  steepest  part  of  the  potential  gradient  Is  between  35°N  and 
55*N  and  thus  the  largest  positive  and  negative  Y-components  of  the 

lunar  tidal  force  exist  In  this  region. 

i 

i 

i 

j  2.5.2  Fx  component. 

The  Fx  component  of  the  lunar  tidal  force  Is  totally  dependent 
on  Terms  2  and  3.  There  Is  absolutely  no  contribution  from  Term  1. 

|  As  discussed  earlier  the  timing  of  the  maximum/minimum  value  of 

potential  at  a  given  longitude  and  latitude  Is  dependent  upon  oq 
(for  Terms  2  and  3).  At  a  given  longitude  where  a  relative  maximum  In 

|  the  east/west  direction  occurs,  the  X-component  would  be  positive  to 

the  east  of  this  longitude  and  negative  to  the  west.  The  opposite 
would  be  true  for  a  longitude  where  a  relative  minimum  occurs. 

Figure  2.7  shows  the  average  Fx  component  of  the  lunar  force  for 
August  1947.  The  potential  Is  a  relative  minimum  approximately  along 
230*W  and  50*W,  thus,  there  Is  a  net  acceleration  towards  these 

i  • 

|  longitudes  In  the  X  direction. 

Table  2.4  shows  a  time  series  (1940-1960)  of  average  monthly 
(August)  values  of  lunar  potential  for  Term  2  plus  Term  3  at  25 *N  for 

|  three  longitudes,  30*E,  75*E,  and  120*E.  Note  that  In  certain  years 

(e.g.,  1947  and  1949)  the  northern  India  location  (25*N,  75*E)  Is  a 


ly  average  of  the  x-component 
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TABLE  2.4 


Average  August  Values  of  Lunar  Potential  (Term  2  +  Term  3)  at  25*N 
for  the  Longitudes  and  Tears  Indicated 
Unite  m2«-2  y  iry-4 


Year 

30  *E 

75*E 

120*E 

1940 

-49.6 

-40.3 

26.6 

1941 

-8.3 

44.6 

48.0 

1942 

41.2 

25.9 

-42.9 

1943 

-16.2 

-68.1 

-49.9 

1944 

-46.7 

-15.8 

42.6 

1945 

9.3 

54.0 

39.8 

1946 

50.3 

5.0 

-75.9 

1947 

-34.7 

-68.9 

-27.4 

1948 

-30.2 

15.0 

53.3 

1949 

30.6 

61.4 

21.5 

1950 

45.7 

-22.9 

-100.0 

1951 

-55.6 

-76.0 

-18.2 

1952 

-40.2 

24.1 

55.0 

1953 

22.1 

59.9 

19.7 

1954 

32.7 

-36.4 

-96.7 

1955 

-57.5 

-59.8 

1.2 

1956 

-22.6 

38.3 

49.9 

1957 

34.5 

50.6 

-11.4 

1958 

13.6 

-49.9 

-76.8 

1959 

-40.9 

-28.6 

22.6 

1960 

-5.1 

50.6 

39.8 

relative  minimum  or  maximum  producing  an  acceleration  either  toward  or 
away  from  the  region. 

2.6  Summary. 

The  nodal  frequency  (.053725  year"*)  dominates  the  spectrum  of 
average  monthly  lunar  potential.  A  frequency  around  .263  year"*  Is 
the  only  other  one  that  explains  more  than  one  percent  of  the 
variance.  Both  of  these  frequencies  are  directly  related  to  some 
aspect  of  taking  the  monthly  average  of  the  fortnightly  tide.  Fig. 

2.8  shows  the  approximate  latitude  of  the  sublunar  point  versus  the 
day  of  the  month  with  rf|nax  equal  to  0.0  or  27.32  and  «max  equal  to 
28.6*M.  The  nodal  cycle  Is  dependent  on  the  range  of  latitude  while 
the  other  lunar  tidal  potential  frequency,  .263  year"*,  is  dependent 
upon  what  the  sublunar  latitude  Is  at  the  beginning  of  the  month.  The 
time  of  maximum  northward  declination  (4mx)  for  the  month  was  the 
parameter  used  to  prescribe  this  latitude  In  Equation  2.28. 

Based  on  the  analysis  In  this  chapter  one  would  expect  that  If 
lunar  tides  In  the  atmosphere  have  an  Impact  on  climatic  Interannual 
variability,  the  two  frequencies  found  to  be  the  most  Important  in  the 
Interannual  variability  of  the  lunar  potential  would  exist  In  the  time 
series  of  climatic  variables.  The  working  hypothesis  then  is  that 
atmospheric  lunar  tides  do  produce  a  significant  signal  In  climatic 
variability.  It  now  remains  to  test  this  hypothesis.  The  next 
chapter  will  be  devoted  to  this  effort. 


LATITUDE 


SUBLUNAR  LATITUDE  VS.  DAY  OF  MONTH 


DAY  OF  MONTH 
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CHAPTER  3 

Periodic  Component:  Test  of  Hypothesis  -  Observed  Response 

3.1  Introduction 

Mo  single  theory  of  climatic  variability  can  explain  all  of  the 
observed  Interannual  variance  because  the  climate  Is  controlled  by 
many  factors.  The  last  chapter  presented  one  forcing  mechanism, 
lunar  tides,  that  was  hypothesized  to  account  for  some  fraction  of 
the  climatic  Interannual  variability.  We  saw  that  two  frequencies 
dominate  the  spectrum  of  average  monthly  lunar  tidal  potential.  The 
frequency  that  explained  the  most  variance  (.0537  year-1)  was 
associated  with  the  nodal  cycle.  The  only  other  spectral  peak  of 
any  consequence  was  at  a  frequency  around  .263  year-1  In  June. 

This  frequency  was  associated  with  the  time  of  the  month  that  the 
maximum  north  or  south  declination  occurred  which  essentially 
determined  what  residual  part  of  the  sidereal  month  (27.3216  days) 
did  not  average  out  to  zero  during  the  30  days  of  June. 

In  this  chapter  the  hypothesis  that  lunar  tidal  forces  have  an 
Impact  on  Interannual  variability  will  be  tested.  To  test  the 
hypothesis  we  must  examine  a  region  that  Is  sensitive  to  this 
forcing.  As  Pearson  (1982)  showed,  a  highly  significant  periodic 
signal  (possibly  related  to  the  Chandler  tide)  In  the  December 
temperatures  In  the  upper  Midwest  was  observable  because  this  area 
lies  close  to  the  boundary  between  cold  Arctic  air  and  milder 
Pacific  air  at  this  time  of  year.  Thus,  a  small  shift  from  the 


normal  position  of  the  Arctic  front  can  have  a  large  Influence  on 
the  mean  December  temperatures.  This  periodicity  In  temperature 
results  from  the  Interaction  of  the  forcing  (the  proximal  cause 
being  pressure  anomalies  In  northern  Canada)  and  the  synoptic 
singularity  which  makes  the  upper  Midwest  particularly  sensitive  to 
the  periodic  signature  of  the  forcing.  In  India,  June  Is  marked  by 
another  climatic  singularity,  the  onset  of  the  monsoon.  As 
discussed  In  Chapter  1,  June  precipitation  In  northern  India  Is 
particularly  sensitive  to  the  timing  of  the  onset  just  as  December 
temperatures  In  the  upper  Midwest  are  sensitive  to  the  timing  of  the 
transition  from  a  late  fall  to  winter  flow  pattern. 

Lunar  tidal  forcing  Is  always  acting  on  the  atmosphere; 
however.  It  Is  likely  that  Its  effect  would  be  most  readily 
observable  in  a  region  that  Is  especially  sensitive  to  a  synoptic 
singularity  like  the  monsoon  front  where  small  departures  from  the 
front's  normal  position  can  have  a  large  Impact  on  the  percent  of 
normal  rainfall  received.  Therefore,  If  lunar  forcing  does  have  an 
effect  on  the  climate,  then  In  this  sensitive  region  one  or  both  of 
the  two  lunar  frequencies  described  In  Chapter  2  might  be  expected 
to  exist  In  the  spectrum  of  Interannual  variability. 

The  lunar  forcing  hypothesis  will  be  tested  by  examining  how 
average  monthly  June  station  pressure,  cube  rooted  precipitation, 
and  300  mb  height  data  respond  at  these  two  hypothesized 
frequencies.  First  the  responses  of  these  three  climatic  variables 
at  the  hypothesized  frequencies  will  be  calculated  to  see  If  the 


signal  Is  strong  enough  to  produce  a  consistent,  synoptic 
climatology  between  them.  If  so,  then  It  remains  to  be  seen  If 
these  three  variables  are  dynamically  related  to  and  In  phase  with 
the  hypothesized  forcing. 

It  should  be  emphasized  that  this  research  Is  dealing  with 
Interannual  variability.  In  northern  India  the  Interannual 
variability  of  June  precipitation  Is  about  38  percent  of  the 
Intra-annual  variance  of  monthly  precipitation  and,  thus,  equals 
only  one  to  two  percent  of  the  daily  variance  (If  one  considers  that 
the  variance  of  monthly  precipitation  amounts  is  about  1/30  that  of 
the  dally  variance  using  the  sampling  distribution  of  means). 
Therefore,  even  If  lunar  tides  account  for  a  very  small  fraction  of 
the  total  variance,  this  forcing  mechanism  becomes  a  candidate 
causal  factor  to  consider  when  trying  to  understand  Interannual 
vari  ability* 

3.1.1  Atmospheric  Susceptibility. 

Under  the  hypothesis  of  an  external  forcing  acting  on  the 
atmosphere,  what  Is  the  nature  of  the  response?  The  relationship 
between  forcing  and  response  can  be  written  simply  as  : 

R  «  kFL(t)  3.1 

where  R  Is  the  atmospheric  response  to  the  forcing, 

k  Is  the  susceptibility  of  the  atmosphere  to  the  forcing, 
and 

Fl  Is  the  lunar  forcing  which  Is  a  function  of  time. 


In  the  discussion  that  follows  the  nomenclature,  FQ5  and  F2g, 
will  be  used  to  represent  the  lunar  forcing  at  the  frequencies  .0537 
year"1  and  .263  year"1,  respectively. 

In  the  case  of  December  temperatures,  described  by  Pearson, 
Equation  3.1  could  be  written  as 

Rtemp  a  KosiZift  *  4)  3.2 

where  Rtemp  temperature  to  the  forcing, 

A  Is  the  amplitude  of  the  response,  and 
cos(2»ft  +  d)  Is  the  periodic  forcing. 

The  amplitude.  A,  of  the  response  Is  the  factor  k  of  Equation  3.1 
and  Indicates  how  susceptible  the  temperature  Is  to  the  forcing. 

Let  us  now  examine  what  kind  of  susceptibility  factor,  k,  might 
relate  Indian  monsoon  climatic  variables  to  the  hypothesized  lunar 
forcing.  The  factor,  k.  In  the  December  temperature  variability  was 
taken  as  a  constant  and  made  the  relationship  between  forcing  and 
response  linear.  However,  considering  the  nonlinearity  of  the 
atmosphere,  this  factor  certainly  does  not  have  to  be  linear.  Let 
us  assume  that  the  atmosphere  Is  more  susceptible  to  lunar  forcing 
at  certain  times  than  at  others.  If  this  Is  so,  then  the 
susceptibility  factor  Is  also  a  function  of  time  and  the  problem  Is 
no  longer  linear,  for  then  we  have: 

R  -  k(t)FL(t).  3.3 

It  Is  reasonable  to  think  of  k  as  being  related  to  the  strength  of 
the  monsoon  circulation,  for  when  the  monsoon  circulation  Is 
particularly  strong,  one  would  expect  the  effects  of  lunar  forcing 


to  be  minimal  because  northern  India  would  no  longer  be  sensitive  to 
small  fluctuations  In  the  timing  of  monsoon  onset  since  with  a 
strong  circulation  the  monsoon  would  arrive  well  before  the  end  of 
the  month. 

Hemispheric  temperature  Is  one  factor  which  has  an  Important 
effect  on  the  Intensity  of  the  monsoon  circulation.  Bryson  and 
Swain  (1981)  compared  a  reconstructed  history  of  monsoon  rainfall  in 
Rajasthan,  India  to  what  was  happening  to  the  temperature  in  the 
high  latitudes  and  found  that  the  period  of  very  heavy  monsoon  rains 
(after  10,800  BP  up  to  3,700  BP)  corresponded  to  a  period  when  the 
region  of  cold  Arctic  climate  (defined  as  the  area  to  the  north  of 
the  tundra/boreal  forest  boundary)  had  retreated  280  km  to  the  north 
of  Its  present  day  position.  During  the  monsoon  drought  period  In 
Rajasthan  from  3,700  BP  to  2,000  BP  the  tundra  region  expanded  to 
the  south  of  Its  present  position.  Other  fluctuations  In  the 
monsoon  rainfall  since  2,000  BP  can  similarly  be  related  to  the 
Arctic  temperatures.  Studies  (e.g..  Van  Loon  and  Williams,  1976; 
Brlnkmann,  1979;  and  Groveman  and  landsberg,  1979)  have  shown  that 
the  temperature  trend  that  occurs  In  the  Arctic  and  Subarctic  region 
determines  to  a  very  large  extent  the  overall  trend  In  Northern 
Hemispheric  temperature.  Thus,  a  definite  relationship  should  exist 
between  Northern  Hemisphere  temperature  changes  and  the  Intensity  of 
the  monsoon  circulation.  Providing  further  support  to  this  Idea 
with  an  emphasis  on  seasonality  Is  a  study  (Kutzbach  and 
Otto-Bllesner,  1982)  which  showed,  by  using  a  low  resolution  general 
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circulation  model,  that  9000  years  ago  when  there  was  estimated  to 
be  7  percent  more  solar  radiation  received  during  the  Northern 
Hemisphere  summer,  the  monsoon  circulation  was  more  Intense.  Thus, 
heml spherically  warm  periods  produce  strong  monsoons;  whereas,  cold 
periods  have  more  winter-1  Ike  circulation  patterns  with  weak  monsoon 
circulations.  Therefore,  If  we  assume  that  the  susceptibility 
factor  Is  dependent  on  the  strength  of  the  monsoon  circulation,  then 
making  It  a  function  of  Northern  Hemisphere  temperature  seems  like 
an  appropriate  step. 

One  source  of  data  to  determine  a  susceptibility  factor  based  on 
Northern  Hemispheric  temperature  comes  from  Jones  et.  al.  (1982). 
They  have  published  a  100  year  time  series  of  mean  monthly  and 
annual  Northern  Hemispheric  temperatures.  While  there  are 
difficulties  with  any  such  data  set,  the  long  term  characteristics 
of  Northern  Hemispheric  temperature  may  be  determined.  Fig  3.1 
shows  the  mean  annual  Northern  Hemslphere  surface  temperature 
anomalies  from  Jones  et.  al.  (op.  clt.).  A  definite  warming  trend 
Into  the  1940 's  and  then  cooling  Into  the  1970's  Is  observed.  As  a 
first  approximation  to  this  temperature  data,  the  best  fitting 
cosine  curve  was  found  to  be: 

TJWK  *  0.37cos[2»( .0057923)t  -  2.464],  3.4 

This  curve  has  a  period  of  about  172.6  years  with  a  maximum  value 
occurring  In  1948.7,  and  It  explains  58  percent  of  the  variance  of 
the  mean  annual  Northern  Hemslphere  surface  temperature  time 
series.  Looking  at  the  time  series  of  average  temperature  for  the 
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Figure  3.1.  The  solid  line  Is  the  Northern  Hemisphere  annual  mean 
surface  temperature  anomalies  from  the  1946-1960  mean 
(after  Jones  et.  al.,  1982).  The  dashed  line  is  the 
best  fitting  cosine  (explaining  58  percent  of  the 
variance)  used  In  the  susceptibility  factor. 


months  April  through  July  (l.e.,  around  the  time  of  the  monsoon 
onset)  the  best  fitting  cosine  Is  very  similar  to  the  one  for  the 
annual  data. 

A  simple  a  priori  hypothesis  on  how  to  Incorporate  this 
Information  Into  a  susceptibility  factor  Is  to  assume  that  when  the 
hemisphere  Is  warmest  the  tidal  forcing  has  no  effect  on  the  monsoon 
circulation.  This  produces  the  following  model  for  the  theoretical 
response, 

R  -  A[Q(t)  -  B]Fl  3.5 

where  B  Is  set  equal  to  1.0,  Q(t)  is  the  cosine  approximation  to  the 
Northern  Hemisphere  surface  temperature  In  Equation  3.4,  and  F|_  Is 
the  lunar  forcing  cosine.  Thus,  the  factor  In  brackets  ranges  from 
a  zero  (wannest  temperatures)  to  minus  two  (coldest  temperatures). 
The  parameter  A  allows  for  the  sign  and  magnitude  of  the  response. 

It  Is  Important  to  emphasize  that  this  form  of  the  susceptibility 
factor  Is  only  a  first  attempt  to  define  this  aspect  of  the 
variability.  It  may  be,  for  example,  that  Q(t)  and/or  the  entire 
susceptibility  factor  should  be  raised  to  some  power.  Further 
research  will  be  required  to  gain  greater  insight  Into  the 
characteristics  of  this  nonlinear  factor. 

Now  the  nature  of  the  theoretical  response  will  be  examined.  As 
Equation  3.5  Indicates,  there  are  three  main  components  to  this 
response: 

1.  The  lunar  forcing  cosine,  FL.  The  frequencies  used 
(.0537  year-*  and  .263  year"*)  and  the  phase  of  each  forcing 


frequency  were  determined  from  the  perlodogram  calculated  In  Chapter 
2  and  shown  In  Figure  3.2.1; 

2.  The  cosine  approximation,  Q(t),  to  the  Northern 
Hemisphere  surface  temperature  which  was  discussed  above;  and 

3.  The  term  B  which  will  be  called  the  susceptibility 
threshold.  The  value  of  B  will  be  empirically  determined  from  the 
data. 

As  noted  above.  If  B  Is  one  then  at  the  warmest  times  the 
response  to  the  forcing  Is  zero.  If  B  Is  less  than  one  then  for  a 
certain  period  of  time  the  response  actually  reverses  sign.  If  B  Is 
greater  than  one  then  the  response  Is  always  susceptible  to  the 
forcing  but  at  varying  degrees  and  It  always  has  the  same  sign.  As 
B  Increases  the  hemispheric  temperature  effect  becomes  less  and  less 
Important.  Thus,  If  B  Is  large  enough,  this  effect  Is  essentially 
put  to  zero  and  only  the  forcing  Itself  Is  Important.  Figure  3.3 
shows  the  effect  that  the  susceptibility  threshold  has.  For  the 
purposes  of  demonstration,  a  forcing  frequency  of  0.20  year"1  and 
a  modulating  period  of  180  years  were  used.  In  all  four  cases  In 
Figure  3.3  year  20  has  a  relative  minimum  In  the  cycle.  However, 
depending  on  the  value  of  B,  year  90  Is  either  In  phase,  zero,  or 
180  degrees  out  of  phase  with  year  20.  Thus,  the  empirically 
determined  value  of  the  susceptibility  threshold,  B,  determines  how 
the  temperature  cosine  (Q(t)  In  Eqn.  3.5)  actually  modulates  the 
lunar  forcing.  We  will  see  that  the  value  of  the  susceptibility 
threshold  has  a  regional  character. 


59 


FIGURE  3.2.1  -  LUNAR  POTENTIAL  (TERM  1) 


FIGURE  3.2.2  -  TMEDLF  RESPONSE 


Figure  3.2.  Perlodograms  of  the  lunar  tidal  potential  (from  Term  1) 
and  the  associated  TMEDLF  response  for  the  years  1899 
through  1960.  The  TMEDLF  response  time  series  was 
calculated  using  the  amplitudes  and  phases  for  the 
frequencies  .0537  year-’  and  .263  year-'  from  the  lunar 
tidal  potential  time  series  (Term  1)  along  with  the 
representative  values  of  the  susceptibility  threshold 
as  shown  In  Equation  3.6. 
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Figure  3.3  Time  series  for  a  response  to  a  forcing  with  a  five  year 
period  and  a  180  year  modulation.  The  only  difference 
between  each  time  series  Is  the  value  of  the  susceptibility 
threshold.  These  plots  show  how  sensitive  the  response  is 
to  the  value  of  the  susceptibility  threshold  in  Eqn  3.5. 
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Since  In  this  model  of  the  response  the  amplitude  of  the  lunar 
forcing  Is  being  modulated  by  a  cosine  derived  from  temperature  data 
and  since  the  characteristics  of  this  modulation  are  governed  by  the 
empirically  determined  susceptibility  threshold  B,  the  expression 
"temperature  modulation,  empirically  determined,  of  the  lunar 
forcing"  (TMEDLF)  Is  appropriate.  Therefore,  throughout  the  rest  of 
this  thesis  the  theoretical  response  will  be  referred  to  as  the 
TMEDLF  response.  It  Is  the  expected  response  based  on  the  above 
theoretical  discussion  of  lunar  forcing  and  atmospheric 
susceptibility  to  this  forcing.  The  TMEDLF  response  will  be 
compared  to  the  actual  response  observed  In  the  pressure, 
precipitation,  and  300  mb  data  to  see  If  this  observed  response  Is, 
In  fact,  similar  to  the  one  predicted  (l.e.,  the  TMEDLF  response). 

A  perlodogram  of  the  forcing  time  series  developed  In  Chapter  2 
(Fig.  3.2.1)  shows  the  large  peak  at  .0537  year-1,  and  the  much 
smaller  one  at  .263  year”*.  We  saw  In  Table  2.2  that  these  two 
frequencies  alone  explained  almost  100  percent  of  the  lunar 
potential  variance.  Therefore,  to  determine  the  nature  of  the 
TMEDLF  response  to  this  forcing,  we  can  multiply  each  periodic 
component  by  the  appropriate  susceptibility  factor  as  follows: 

R  -  A1CQ(t)  -  B1]cos(2»f1t  ♦  dx)  ♦ 

Ag[Q(t)  -  B2]cos(2irf2t  ♦  rf2)  ♦  C  3.6 
where  Q(t)  Is  the  cosine  approximation  to  the  Northern 
Hemispheric  surface  temperature. 


fj  Is  the  frequency,  .0537  year-*,  in  the  lunar 
tidal  potential, 

A^,  are  the  amplitude  and  phase  of  f^, 

Bj  is  the  susceptibility  threshold  for  the  response 
at  the  f^  frequency, 

f2  is  the  frequency,  .263  year-1,  in  the  lunar 
tidal  potential, 

A2,  d2»  B2  are  t0  f2  as  Ai»  ^1*  Bi  are  t0 
fj,  and 

C  Is  the  mean  potential  value. 

There  Is  no  reason  to  assume  that  the  susceptibility  factor  Is 
the  same  for  both  frequencies  and.  In  fact,  we  will  see  that  the 
empirically  determined  value  for  term  B  Is  higher  for  the  response 
at  .263  year"1  than  It  Is  for  the  response  at  .0537  year-1.  If 
we  set  B^  equal  to  0.93  and  B2  equal  to  1.8,  which  are  values 
found  In  the  regressions  of  the  climatic  variables  (see  Section 
3.5),  then  the  periodogram  of  Fig.  3.2.2  results  with  the  peaks  at 
.0484  year-1  and  .263  year-1  explaining  32.2  and  20.6  percent  of 
the  TMEDLF  response  variance,  respectively.  There  are  several 
Important  points  to  note: 

a.  The  peak  at  the  nodal  cycle  no  longer  exists  but  has  shifted 
to  the  lower  frequency  of  .0484  year-1  and  Its  amplitude  has 
decreased  with  an  Indication  of  band  splitting.  This  Is  the  result 
of  the  amplitude  modulation. 
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b.  When  a  representative  value  for  B  Is  used  for  both  lunar 
frequencies,  even  though  In  the  lunar  potential  time  series  the 
amplitude  of  .263  year"1  Is  considerably  less  than  the  amplitude 
of  the  nodal  (see  Fig.  3.2.1),  the  spectral  characteristics  of  the 
TMEDLF  response  no  longer  reflect  this  difference  but  Instead  the 
amplitudes  of  the  two  frequencies  are  fairly  close. 

c.  The  peak  at  .263  year"1  Is  not  shifted.  It  exists  In  both 
the  forcing  and  the  TMEDLF  response. 

3.1.2  Methodology  used  to  test  the  hypothesis. 

In  order  to  determine  the  behavior  of  the  station  pressure, 
precipitation  and  upper  air  data  using  the  concept  of  atmospheric 
susceptibility  to  external  forcing,  the  MACC  nonlinear  regression 
program  NREG  was  used.  The  regression  equation  was: 

P  .  A[cos(2*fjt  ♦  dx)  -  B]cos(2irf2t  +  rf2)  ♦  C  3.7 
where  P  Is  the  predicted  value  of  the  climatic  element  (e.g., 
station  pressure)  from  the  regression, 

A  provides  a  magnitude  and  sign  for  the  response, 
f^,  dj  are  the  fixed  frequency  and  phase  of  the  best 
fitting  cosine  to  the  Northern  Hemisphere  temperature 
time  series, 

B  Is  the  susceptibility  threshold  (see  below  for 
discussion), 

f2  Is  the  forcing  frequency  (.053725  year"1  or  .263 

-1\ 
year  x), 
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d2  Is  the  phase  of  the  forcing  and  Is  allowed  to  float, 
and 

C  Is  the  mean. 

The  NREG  routine  returns  the  values  for  the  four  parameters 
(A,B,C,  and  d2)  plus  the  predicted  values  for  the  years  of  the 
regression.  The  susceptibility  threshold  B  Is  not  fixed  In  the 
regression  so  that  we  may  see  how  sensitive  the  particular  climatic 
variable  Is  to  the  effects  of  hemispheric  warming  and  cooling.  The 
phase  of  the  forcing,  Is  allowed  to  float  to  test  whether  or 
not  the  pressure,  precipitation  and  upper  air  data  are.  In  fact.  In 
phase  with  the  hypothesized  lunar  forcing. 

To  obtain  the  percent  explained  variance,  the  square  of  the 

correlation  coefficient  (R2)  between  the  observed  and  the 

predicted  values  was  calculated.  In  the  analysis  that  follows  the 
2 

maps  of  R  will  be  presented  along  with  maps  of  the  predicted 
values  for  certain  years.  In  some  cases  the  map  of  the 
susceptibility  threshold  will  also  be  shown.  The  factor  A  will  not 
be  shown  since  Its  magnitude  does  not  necessarily  Indicate  the 
magnitude  of  the  response  to  the  forcing  since  It  must  also  be 
multiplied  hy  the  term,  B.  Thus,  A  means  different  things  at 
different  stations  and  Is  not  easily  comparable.  The  phase  maps  of 
d2  will  also  not  be  shown  since  the  term  A  can  have  a  negative 
sign,  so  at  the  stations  where  this  occurs  the  phase  Is  180*  from 
what  It  should  be  to  compare  It  to  the  other  stations.  Instead  of 
adding  180*  to  the  phase  In  those  cases  and  then  plotting  the  data, 
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•the  phasing  Information  will  be  derived  from  looking  at  the  actual 
anomaly  maps  for  the  Important  years  In  the  cycle. 

In  order  to  test  the  significance  of  the  regression,  an  F  ratio 
test  can  be  applied  for  a  3  parameter  model  as  follows: 

F  «  R2/3  3.8 

(l-R2)/N-4 

where  R  Is  the  correlation  between  the  predicted  values  from  the 

NREG  regression  and  the  actual  observed  values,  and  N  Is  the  number 

of  data  points.  Testing  the  significance  at  the  95  percent  level 

2 

establishes  the  F-ratlo  for  the  given  degrees  of  freedom.  The  R 

value  can  then  be  written  as: 

R2  -  D/U  +  D)  3.9 

2 

where  D  *  3F/(N-4).  Thus,  apy  R  value  greater  than  11.8  percent 

Is  significant  at  the  95  percent  level  for  66  years  of  data.  For  27 
2 

years  of  data  R  must  be  greater  than  26.5  percent. 

In  the  discussion  that  follows  the  phrases,  "susceptibility 
factor",  "susceptibility  parameter",  and  "susceptibility  threshold" 
will  be  used.  The  first  pertains  to  the  entire  factor, 

A[cos(2vfjt  ♦  6^)  -  B],  the  second  to  the  cosine  term  of  this 
factor  and  the  latter  to  the  term,  B.  The  phrase  "TMEDLF  response" 
will  mean  the  forcing  frequency  multiplied  by  the  susceptibility 
factor.  Also,  the  nomenclature  kFQg  and  kFgg  will  refer  to  the 
nonlinear  model  of  Equation  3.7  depending  on  whether  f ^  In  this 
equation  Is  .0537  year"*  or  .263  year"*,  respectively. 

It  will  be  necessary  to  correlate  the  responses  observed  In 
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pressure,  precipitation  and  the  upper  air  data  with  the  tidal 
forcing  to  see  If  they  are,  In  fact,  in  phase.  The  correlation 
between  two  time  series  generated  by  cosines  with  Identical 
frequencies  but  different  phases  Is  merely  the  cosine  of  the  phase 
difference.  However,  our  situation  Is  slightly  more  complicated 
since  the  forcing  cosine  is  multiplied  by  the  susceptibility 
factor.  Therefore,  a  time  series  of  the  TMEDLF  response  (l.e.,  the 
susceptibility  factor  multiplied  by  the  lunar  forcing)  must  be 
calculated  and  correlated  with  the  P  time  series  of  predicted  values 
from  the  nonlinear  regression  from  Equation  3.7.  The  susceptibility 
threshold  used  In  the  TMEDLF  response  (Eqn.  3.5)  will  be  the 
empirically  determined  value  for  the  station  under  consideration. 

At  times  the  TMEDLF  response  will  be  correlated  with  the  average 
anomaly  for  a  group  of  stations.  When  this  occurs  the 
susceptibility  threshold  used  In  the  TMEDLF  response  will  be  the 
averse  value  of  the  susceptibility  threshold  for  the  stations  being 
considered.  !f  a  station  has  a  susceptibility  threshold  that  Is 
greater  than  5.0,  however,  then  that  term  will  not  be  Included  In 
the  threshold  average,  but  the  station's  anomaly  will  still  be 
used.  The  reason  for  this  Is  that  a  threshold  value  that  high 
Indicates  the  station's  response  Is  not  dependent  upon  the 
temperature  modulation  of  the  forcing.  The  phasing  Information  for 
the  forcing  was  obtained  from  the  lunar  tidal  potential  time  series 
for  25*N,  75*E  (a  location  In  northern  India).  However,  as  we  saw 
In  Chapter  2,  there  Is  only  a  latitudinal  dependence  for  the  two 


frequencies  being  considered  so  the  longitude  of  the  location  Is  not 
Important.  The  Fqj  cycle  has  a  maximum  In  June  1913,  for  the 
Moon's  ascending  node  was  at  the  verna'  equinox  at  that  time.  The 
phase  for  the  F2g  cycle,  taken  from  the  perlodogram  of  Term  1,  has 
a  maximum  In  1900.65.  For  the  sake  of  consistency  In  these 
correlation  calculations,  the  sign  of  the  TMEDLF  response 
coefficient.  A,  In  Eqn.  3.5  was  taken  to  be  that  from  the  regression 
of  the  300  mb  data  at  25*M,  75*E,  for  as  we  shall  see  the  upper  air 
anomalies  over  northern  India  are  an  Important  factor  in  the 
synoptic  climatology  of  the  monsoon  onset. 

3.2  Data 

3.2.1  Station  Pressure. 

I  have  used  WMO  station  pressure  data  at  40  stations  In  India 
and  Pakistan  for  this  study.  Figure  3.4  shows  the  locations  and 
Table  3.1  provides  the  WMO  numbers  and  names  of  these  stations. 
Station  pressure  was  used  rather  than  sea  level  pressure  for  two 
main  reasons: 

1.  The  data  record  length  for  most  stations  was  longer  for 
station  pressure  than  for  sea  level  pressure. 

2.  Station  pressure  Is  the  actual  observation  and  has  not 
been  subjected  to  the  reduction  to  sea  level  formula.  Since  we  are 
looking  at  periodicities  In  the  data  It  seemed  best  to  use  unaltered 
data*  The  resulting  anomalies  are  departures  from  the  station  mean 
and,  therefore,  even  If  the  means  are  quite  different  due  to  station 
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TABLE  3.1 

WMO  Numbers  and  Names  of  Stations  In  Figure  3.4  (Station  Pressure) 


4  number 

Name 

WMO  Number 

1 

Lahore  City 

416400 

2 

Quetta 

416610 

3 

Karachi 

417820 

4 

Simla 

420830 

5 

Ludhiana 

420990 

6 

Mukteswar 

421470 

7 

B1 kaner 

421650 

8 

Agra 

422610 

9 

Darjeeling 

422950 

10 

Jodhpur 

423390 

11 

Jaipur 

423480 

12 

Darbhanga 

423910 

13 

Dhubrl 

424040 

14 

Kota 

424510 

15 

A1 1 ahabad 

424750 

16 

Daltonganj 

425870 

17 

Dumka 

425990 

18 

SI 1 char 

426190 

19 

Sagar 

426710 

20 

Dwarka 

427310 

21 

Indore 

427540 

22 

Calcutta 

428070 

23 

Nagpur 

428670 

24 

Veraval 

429090 

25 

Akola 

429330 

26 

Cuttack 

429700 

27 

Bombay 

430570 

28 

Poona 

430630 

29 

Begampet 

431280 

30 

Vlshakapatnam 

431490 

31 

Madras 

432790 

32 

Bangalore 

432950 

33 

Amlnl 

433110 

34 

Port  Blair 

433330 

35 

Kodlakanal 

433390 

36 

Fort  Cochin 

433510 

37 

Pambam 

433630 

38 

Trivandrum 

433710 

39 

Srinagar 

435400 

40 

Leh  Kashmir 

435440 

*  *  » *  • 
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elevation  differences,  the  anomalies  are  comparable. 

There  Is  a  great  deal  of  missing  data  after  1960  In  most  of  the 
records.  Therefore,  the  time  period  that  will  be  used  covers 
1895-1960. 

3.2.2  Precipitation 

As  with  station  pressure  the  data  comes  from  the  UMO  and  the 
locations  of  the  39  stations  used  are  shown  In  Figure  3.5.  Table 

3.2  provides  the  WMO  numbers  and  names  of  these  stations. 

The  data  record  for  these  stations  goes  through  1982;  however, 
to  be  comparable  to  the  pressure  data  set,  the  years  1895  -  1960 
will  be  used  In  the  regressions. 

The  frequency  distribution  of  rainfall  amounts  for  a  given  month 
can  be  highly  skewed  especially  In  the  drier  regions  of  India.  If 
the  data  are  not  normally  distributed,  standard  significance  tests 
do  not  apply.  To  reduce  the  skewness,  the  precipitation  data  were 
subjected  to  a  cuberoot  transform.  Thl s  makes  the  frequency 
distribution  at  a  given  station  more  nearly  normal. 

3.2.3  300  mb  Heights. 

300  mb  heights  are  available  on  tape  every  five  degrees  of 
latitude  (90*N  to  20*N)  and  longitude  for  the  years  1950  through 
1977.  The  data  for  the  year  1958  are  missing  so  there  are  only  27 
years  of  data. 
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TABLE  3.2 

WMO  Numbers  and  Names  of  Stations  In  Figure  3.5  (Precipitation) 


4  number 

Name 

WMO  Number 

1 

Peshawar 

415300 

2 

Lahore  City 

416400 

3 

Quetta 

416610 

4 

Kalat 

416960 

5 

Hyderabad 

417650 

6 

Karachi 

417820 

7 

Mukteswar 

421470 

8 

Bikaner 

421650 

9 

New  Delhi 

421820 

10 

Agra 

422610 

11 

Jodhpur 

423390 

12 

Darbhanga 

423910 

13 

Dhubrl 

424040 

14 

Allahabad 

424750 

15 

Dumka 

425990 

16 

Ahmadabad 

426470 

17 

Sagar 

426710 

18 

Dwarka 

427310 

19 

Indore 

427540 

20 

Calcutta 

428070 

21 

Nagpur 

428670 

22 

Veraval 

429090 

23 

Akola 

429330 

24 

Jagdalpur 

430410 

25 

Bombay 

4X570 

26 

Poona 

430630 

27 

Vlshakapatnam 

431490 

28 

Masullpatam 

431850 

29 

Bel gaum 

431970 

X 

Madras 

432790 

31 

Mangalore 

432830 

32 

Bangalore 

432950 

33 

Port  Blair 

433330 

34 

Pamban 

43X30 

35 

Mfnfcoy 

4X690 

36 

Trivandrum 

433710 

37 

Trine  omalee 

434180 

38 

Colombo 

434660 

39 

Sri  nag ar 

435400 
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3.3  Response  to  lunar  forcing  at  .0537  year"1. 

3.3.1  Station  Pressure. 

In  this  section  I  will  present  an  analysis  of  the  June  monsoon 

circulation  as  exhibited  by  the  surface  pressure  field  by  examining 

the  June  anomaly  patterns  for  each  year  produced  by  the  hypothesized 

model  of  Equation  3.7., 

2 

The  R  map  (Fig  3.6)  shows  two  major  regions  of  Interest  -  one 
In  northeastern  India  and  another  In  western  India.  The  phase  map 
for  the  nodal  cycle  (not  shown)  Indicates  that  the  pressure 
variation  In  northeastern  India  Is  close  to  180*  out  of  phase  with 
the  region  In  western  India. 

Looking  at  the  two  areas  outlined  in  Figure  3.7,  we  see  that  the 
susceptibility  threshold  Is  greater  than  one  for  most  of  the 
stations  of  Interest  In  the  west  while  It  Is  less  than  one  In  the 
northeast.  No  attempt  has  been  made  to  analyze  these  data,  since 
the  values  vary  so  much.  In  the  northeast,  then,  as  the  hemisphere 
warms,  the  sign  of  the  response  anomaly  reverses.  Thus,  for  a  time 
the  west  and  northeast  anomalies  are  In  phase.  The  Impact  of  this 
will  be  discussed  below. 

The  oscillating  nature  of  the  anomaly  pattern  at  the  modulated 
nodal  frequency  can  be  seen  by  examining  the  years  1895,  1903,  1913, 
and  1922  shown  In  Fig  3.8.  These  years  were  chosen  because  they  are 
the  Junes  that  are  at  the  extremes  of  the  TMEDLF  response  at  the 
nodal  cycle.  There  was  a  maximum  In  the  TMEDLF  response  In  1894; 
however,  since  the  data  set  began  In  1895,  that  year  was  used. 


Figure  3.6.  The  percent  of  th« 
by  the  kFgg  model 
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Figure  3.7.  The  value  of  the  susceptibility  threshold  from  the 
regression  of  the  station  pressure  data  against  the 
kpQg  model  (units  are  value  x  10-2).  The  two  areas 

outlined  show  that  the  threshold  values  are  different 
In  northeastern  India  from  what  they  are  In  western 
India. 


Figure 
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Also,  the  other  years  are  not  necessarily  a  maximum  or  minimum  In 
the  nodal  cycle  but  may  be  a  year  off  because  of  the  effect  of  the 
amplitude  modulation.  For  example,  1904  was  an  extreme  In  the  nodal 
cycle  while  1903  was  the  extreme  In  the  TWEDLF  response  (see  Table 
3.1).  In  1895  and  again  In  1913,  a  large  positive  anomaly  dominated 
northeastern  India,  extending  northwestard  along  the  Himalayas, 
whereas  a  large  negative  anomaly  occurred  over  western  India.  In 
the  years  1903  and  1922  the  reverse  pattern  occurred.  We  have 
postulated  that  the  tidal  effects  are  less  Important  In  warmer  years 
and  we  see  that  the  magnitude  of  the  anomalies  decreases  from  1895 
to  1922  -  a  period  of  warming  In  the  cosine  approximation  to  the 
Northern  Hemisphere  temperatures.  The  axis  of  the  reversing  pattern 
lies  just  about  along  the  average  position  of  the  monsoon  trough  In 
June  (see  Figure  3.9)  with  the  anomaly  center  to  the  north  of  this 
axis  lying  along  the  Himalayas.  Research  (e.g.,  Raghavan,  1973)  has 
shown  that  when  the  pressure  anomaly  along  the  Himalayas  Is 
negative,  breaks  occur  In  the  monsoon. 

The  effect  of  a  large  positive  anomaly  on  the  west  coast  would 
be  to  decrease  the  strength  of  the  southwesterly  monsoon  current, 
thereby  decreasing  monsoon  rainfall.  An  Inference  about  the  wind 
flow  can  also  be  made  from  these  unnormal  1 zed  pressure  anomaly  maps 
(Kutzbach,  1970).  The  high  and  low  pressure  anomalies.  In  the  years 
depicted,  create  a  strong  NE/SW  pressure  gradient  across  the 
Gangetlc  Plain  of  India.  Since  we  are  dealing  with  monthly  averages 
here,  geostrophlc  equilibrium  can  be  assumed  (with  some  allowance 


Figure  3.9.  The  average  position  of  the  monsoon  trough  (solid  line) 
for  the  years  1899  through  1977  taken  from  the  grldded 
Northern  Hemisphere  sea  level  pressure  data  tape  which 
has  been  updated  by  Trenberth  and  Paollno  (1980).  The 
pivotal  axis  of  the  nodal  cycle  anomaly  oscillation  Is 
shown  with  a  dashed  line. 


for  boundary  layer  friction).  In  the  years  when  the  high  pressure 
anomaly  Is  located  along  the  west  coast,  the  geostrophlc  wind  is 
directed  from  the  dry  deserts  of  northwest  India  and  Pakistan  onto 
central  India  adding  further  to  the  dryness  that  has  already  been 
discussed  In  conjunction  with  this  pattern. 

We  are  dealing  with  the  anomalies  In  the  single  month  of  June. 
However,  It  Is  Interesting  to  note  that  Bhalme  and  Mooley  (1980)  In 
a  study  of  the  seasonal  surface  pressure  anomalies  during  drought 
and  flood  years,  found  strikingly  similar  patterns  as  shown  In  Fig. 
3.10.  They  were  not  dealing  with  ary  periodic  component  of  the  data 
In  their  analysis  but  with  the  total  variability  for  the  entire 
monsoon  season.  The  magnitude  of  their  west  coast  departures  (-0.6 
mb  to  +0.8  mb)  are  comparable  to  the  magnitude  of  the  anomalies  from 
this  periodic  model  (-1.0  mb  to  +0.8mb).  This  Indicates  that 
certainly  some  of  the  variance  associated  with  drought  and  flood 
years  Is  associated  with  this  periodic  component  of  the  variability. 

In  order  to  establish  the  phasing  of  this  periodic  oscillation, 
let  us  examine  the  average  pressure  anomaly  of  the  following  seven 
western  stations  and  eight  northeastern  stations: 


WEST 

NORTH 

WMO 

Fig.  3.4  no. 

UNO 

Fie 

l.  3.4 

47711  -  Dwarka 

20 

47251 

-  Agra 

8 

42909  -  Yeraval 

24 

42475 

-  Allahabad 

15 

43057  -  Bombay 

27 

42587 

-  Daltonganj 

16 

43063  -  Puna 

28 

42970 

-  Cuttack 

26 

42933  -  Akola 

25 

42807 

-  Calcutta 

22 

42754  -  Indore 

21 

42599 

-  Dumka 

17 

42867  -  Nagpur 

23 

42391 

-  Darbhanga 

12 

42099 

-  Ludhiana 

5 

8 


The  difference  between  the  two  averages  (North  minus  West)  will 
provide  a  measure  of  the  strength  of  the  NE/SW  pressure  gradient 
discussed  above.  The  difference  was  calculated  so  that  It  would  be 
negative  for  dry  years.  Table  3.3  provides  a  comparison  of  this 
difference  for  the  years  surrounding  the  maximum  effect  of  the 
TMEDLF  response  (l.e.,  the  susceptibility  factor  times  the  forcing 
cosine).  One  thing  to  note  Is  that  the  pressure  anomaly  difference 
Is,  In  fact,  strongest  at  the  beginning  of  the  data  record  which  Is 
well  Into  the  coldest  period  of  the  best  fit  Northern  Hemslphere 
cosine  curve.  Thus,  the  actual  atmospheric  response  Is  much 
stronger  when  the  hemisphere  Is  coldest.  Also,  the  year  of  the 
maximum  pressure  anomaly  difference  agrees  with  the  year  of  the 
maximum  TfCDLF  response  and  the  correlation  of  t<ie  66  year  time 
series  of  the  pressure  anomaly  difference  and  the  TMEDLF  response  Is 
-.972  when  using  the  average  susceptibility  threshold  of  .94.  The 
minus  sign  of  the  correlation  Is  a  result  of  the  way  the  difference 
In  anomalies  was  calculated.  Thus,  the  actual  atmospheric  response 
Is  very  much  In  phase  with  the  TMEDLF  response.  Note  also  that  even 
though  1904  was  the  year  that  the  lunar  nodal  cycle  produced  a 
mlnlmun  In  the  lunar  tidal  potential  that  when  the  susceptibility 
factor  Is  applied  to  the  forcing  the  minimum  occurs  In  1903  as  It 
also  does  In  the  actual  response. 


The  fact  that  the  susceptibility  term  Is  less  than  one  In 
northeastern  India  means  that  once  the  temperature  cosine  reaches 
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TABLE  3.3 


Pressure  Anomaly  Difference  Compared  to  the  TMEDLF  Response 


TKEDLF  Response 
Threshold  (B-0.94) 


Year 

Difference  (mb) 

(relative 

1895 

1.38 

-128.9 

1896 

1.16 

-109.8 

1902 

-0.96 

90.1 

1903 

-1.07 

100.7 

1904 

-1.04 

99.4 

1905 

-0.91 

87.1 

1911 

0.57 

-58.2 

1912 

0.66 

-67.0 

1913 

0.67 

-67.5 

1914 

0.61 

-60.4 

1921 

-0.29 

37.1 

1922 

-0.31 

37.7 

1923 

-0.30 

34.0 

1924 

-0.27 

27.0 

Note:  The  pressure  anomaly  difference  Is  calculated  by  talcing  the 
difference  of  the  average  anomaly  for  8  northeastern  stations  and  7 
western  stations,  the  TtEDLF  response  uses  the  average 
susceptibility  threshold  for  the  15  stations. 
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that  value  (on  Its  way  to  a  maximum  of  one  at  the  wannest  part  of 
the  cycle)  the  sign  on  the  nodal  forcing  cosine  reverses.  Thus, 
where  the  signs  of  the  pressure  anomalies  In  the  northeast  and  In 
the  west  were  opposite  In  the  early  part  of  the  record,  they  become 
the  same  as  the  hemisphere  warms.  This  greatly  decreases  the  NE/SW 
pressure  gradient  making  the  effect  of  the  response  at  the  nodal 
forcing  f requency  close  to  zero.  This  effect  can  be  seen  by  looking 
at  the  anomaly  map  for  1941  (Fig.  3.11)  where  the  northeastern 
pressure  anomaly  has  a  maximum  of  +.36  mb  (and  +.56  mb  further  to 
the  northwest)  and  the  western  anomaly  has  a  maximum  of  +.37  mb. 

Calculation  of  the  geostrophlc  wind  for  some  of  these  anomaly 
patterns  shows  the  magnitude  of  the  Interannual  variability  Included 
In  this  nonlinear  periodic  model.  For  example.  In  1895  the  pressure 
difference  between  Jaipur  and  Agra  (station  numbers  11  and  8, 
respectively.  In  Fig.  3.4)  was  0.86  mb.  Between  these  two  stations 
lies  the  pivotal  line  of  the  oscillating  pressure  anomalies.  Using 
the  geostrophlc  wind  equation, 

Vg  >  -(l/pfhp/an,  3.10 

where  f  Is  the  corf  oils  parameter, 
p  Is  the  air  density,  and 
ap/an  Is  the  pressure  gradient, 

we  get  Vg  equal  to  4.99  m/s  from  the  south.  In  1955  the  pressure 
difference  between  Begampet  and  VI shakapatnam  (numbers  29  and  30, 
respectively)  was  0.98  mb  which  yields  a  geostrophlc  wind  of  3.49 
m/s  from  the  south.  Crutcher  (1966)  found  that  the  average 


Meridional  component  In  this  region  Is  zero  for  June  with  a  standard 
deviation  of  about  4  kts  (or  2.1  m/s).  This  Indicates  that  the 
Meridional  component  Is  just  as  likely  to  be  from  the  north  as  from 
the  south  and  that  a  large  part  of  that  variation  can  be  explained 
by  the  model  under  discussion. 

3.3.2  Precipitation. 

Having  examined  the  response  of  station  pressure  using  the 
hypothesized  model  (Eqn.  3.7)  and  having  briefly  discussed  the 
pressure  anomaly  pattern's  Impact  on  the  precipitation,  let  us  look 
at  the  result  of  the  nonlinear  regression  of  precipitation  for  the 
same  period  (1895  -  1960)  against  the  same  model. 

The  percent  variance  map  Is  given  In  Fig  3.12.  The  anomaly  maps 
for  the  years  1895,  1903,  1913,  and  1922  are  shown  In  Fig  3.13.  As 
the  previous  discussion  Indicated,  years  1895  and  1913  should  be  wet 
In  the  northwestern  half  of  India  as  they  are.  At  the  opposite 
extreme  of  the  cycle  the  years  1903  and  1922  are  dry.  The  decrease 
In  the  magnitude  of  the  anomalies  due  to  the  amplitude  modulation  Is 
definitely  noticeable  by  1922. 

Rao  (1976)  has  observed  that,  during  breaks  In  the  monsoon, 
rainfall  still  occurs  In  northeast  India  and  In  southern  peninsular 
India.  As  noted  In  the  last  section,  the  deficient  rainfall  of  1903 
and  1922  Is  associated  with  a  pressure  pattern  observed  during 
breaks  In  the  monsoon  and,  consistent  with  this,  the  precipitation 
anomalies  are  positive  In  northeast  and  southern  India. 


Figure  3.13.  Normalized  anomaly  maps  of  precipitation  for  the  years 
Indicated  from  the  kFn(-  model  (units  are  standard 
deviations  x  1(H). 
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Looking  again  at  Fig.  3.8  for  the  years  1895  and  1903,  we  see 
that  Sagar  (station  19  for  pressure  and  17  for  precipitation)  lies 
In  the  region  where  the  effects  of  the  reversing  flow  pattern  should 
be  felt  the  most;  and  we  find  that  the  correlation  between  Sagar 's 
precipitation  response  and  the  TMEDLF  response  Is  -.998.  This  very 
high  correlation  Is  due  to  the  fact  that  the  phase  of  the  forcing 
cosine  In  the  TMEDLF  response  and  the  phase  of  the  observed  response 
at  Sagar  are  almost  Identical.  Thus,  the  lunar  nodal  cycle  Is 
exactly  In  phase  (disregarding  the  sign  of  the  correlation)  with 
Sagar' s  response. 

In  order  to  gain  a  perspective  on  just  how  Important  the 
response  Is  at  this  frequency  In  terms  of  the  total  variability  of 
June  rainfall  let  us  turn  for  a  moment  to  an  eigenvector  analysis. 
Campbell  et.  al.  (1983),  hereafter  referred  to  as  CBB,  have  shown 
that  the  spectrum  of  the  coefficients  of  the  first  eigenvector  of 
June  cuberooted  precipitation  Is  dominated  by  a  frequency  at  .0482 
year"1  and  another  at  .262  year"1  as  shown  In  Figure  3.14  which 
Is  taken  from  CBB.  This  eigenvector  explains  20.3  percent  of  the 
total  Interannual  variance  of  Indian  precipitation  for  June  with  the 
frequency  .0482  year"1  explaining  17.2  percent  and  the  frequency 
.262  year"1  explaining  15.1  percent  of  the  variance  of  this 
eigenvector's  coefficients.  These  are  the  two  main  frequencies  that 
result  when  the  susceptibility  factor  Is  applied  to  the  lunar 
forcing  as  Figure  3.2.2  shows  Indicating  that  the  lunar  forcing 
frequencies  have  a  significant  role  In  explaining  the  Interannual 
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perlodogram  of  the  coefficients  of  the  first  eigenvector  of  June  cuberooted 
eclpitatlon  (1895-1975).  The  two  dominant  peaks  are  at  .0482  yr-1  and 


variability  of  June  precipitation,  especially  In  northern  India 
where  this  eigenvector  has  the  largest  component  values. 

The  second  eigenvector  ( non-normal  1 zed)  of  station  pressure  Is 
also  dominated  by  a  frequency  around  .0482  year'1,  which  explains 
30.5  percent  of  the  variance.  This  second  eigenvector  explains  17.4 
percent  of  the  total  June  station  pressure  variance.  The 
eigenvector  spatial  patterns  for  both  parameters  are  similar  to  the 
responses  observed  from  the  nonlinear  model.  The  pressure 
eigenvector  component  values  are  opposite  In  sign  In  the  northeast 
from  those  In  the  west  and  the  precipitation  component  values  are 
largest  In  the  northwestern  area  of  India  with  opposite  signs  in  the 
northeast  and  In  southern  peninsular  India.  The  frequencies  and 
their  phases  (from  a  perlodogram)  for  the  two  eigenvector  analyses 
ar*»* 

Frequency  (year”1)  Phase  (degrees)  Period 
Precipitation  .04815  -157.0  1895-1975 

Station  pressure  .04848  -154.7  1895-1960 

The  correlation  (over  the  period  1895-1960)  between  the  two 
resulting  time  series  Is  0.993.  Thus,  at  this  frequency,  which  Is 
predicted  based  on  the  concept  of  atmospheric  susceptibility  to  the 
lunar  forcing,  the  pressure  and  precipitation  responses  are  almost 
perfectly  correlated. 

Note:  All  eigenvectors  used  above  represent  significant 
geophysical  signals  at  the  95  percent  level,  using  the  same 
conservative  significance  test  from  Overland  and  Prelsendorfer 


(1982)  'that  was  used  for  the  eigenvector  analysis  discussed  In  C6B. 

As  we  did  with  pressure,  to  get  a  measure  of  the  relationship 
between  the  forcing  and  the  actual  precipitation  response,  the 
correlation  between  the  THEDLF  response  and  the  average 
precipitation  anomaly  for  12  northwestern  India  stations  was  -.971 
when  using  the  average  susceptibility  threshold  of  .969.  This  is 
not  surprising  considering  how  closely  related  the  preclpltalon  and 
pressure  responses  are. 

Bryson  (1975)  published  a  graph  showing  the  percent  of  stations 
In  northwestern  India  having  less  than  half  of  normal  rainfall  for 
the  period  1900-1972  (Fig  3.15).  The  hypothesized  forcing  model 
(Eqn.  3.7)  explains  25.1  percent  of  the  variance  of  this  graph.  The 
susceptibility  threshold  from  the  regression  Is  0.89.  The 
regression  of  the  graph  against  just  the  Jones  et.  al.  (1982)  annual 
data  explains  31.1  percent  of  the  variance.  Combining  both  of  these 
predictors  Into  a  single  model  (the  dashed  line  In  Figure  3.15) 
explains  51.2  percent  of  the  variance  Indicating  some  shared 
variance  but  not  very  much.  Thus,  northern  Indian  rainfall  Is 
highly  dependent  upon  both  the  Northern  Hemisphere  temperature  and 
the  amplitude  modulated  nodal  cycle. 

As  noted  In  Chapter  2,  Bhalme  and  Mooley  (1981)  studied  cyclic 
variations  In  the  Flood  Area  Index  (FAI)  for  India  (see  Fig.  3.16 
for  the  FAI>1  and  FAI>2  time  series)  and  concluded  that  there  is 
strong  statistical  evidence  for  the  relationship  between  the  double 
sunspot  cycle  and  the  areal  extent  of  flooding  over  India  for  their 


Figure  3.15.  The  solid  line  Is  the  percent  of  the  stations  in 
northern  India  with  less  than  half  of  the  normal 
rainfall  In  a  given  year  (running  ten  year  averages). 

The  dashed  line  represents  the  predicted  values  from  a 
regression  of  the  data  against  both  the  kFgg  model  and 
the  Northern  Hemisphere  annual  mean  surface  temperature. 
The  dashed  curve  explains  over  51  percent  of  the  variance 
The  data  for  the  solid  curve  were  taken  from  Bryson,  1975. 


Figure  3.16.  Flood  Area  Index  time  series  for  the  Indicated  flood 

Intensities  (after  Bhalme  and  Mooley,  1981).  The  dashed 
curves  represent  the  predicted  values  of  a  regression  of 
these  data  against  the  22  year  double  sunspot  cycle, 
while  the  solid  curves  represent  the  predicted  values  of 
a  regression  of  these  data  against  the  kFftC  model. 


spectral  analysis  had  a  peak  at  band  4  (.0455  year-1  which  Is  a  22 
year  period).  A  perlodogram  analysis  showed  that  the  spectral  peak 
was  at  a  slightly  higher  frequency  than  their  harmonic  analysis 
could  resolve.  The  perlodogram  frequency  was  about  .048  year-1 
which  Is  the  same  frequency  found  In  the  linear  eigenvector  analysis 
of  June  station  pressure  and  precipitation. 

When  the  two  FAI  Indices  were  regressed  against  the  kFQ5 
model,  the  FA1>1  regression  explained  17.8  percent  of  the  variance 
compared  to  the  12  percent  from  their  harmonic  analysis.  The 
respective  percents  for  the  FAI>2  regression  were  26.9  and  16.0. 

One  would  expect  the  kFgg  regression  model  to  explain  somewhat 
more  of  the  variance  since  It  has  three  parameters  (which  lead  to 
four  predictors  In  the  linear  equivalent)  whereas  the  harmonic 
analysis  only  has  two  predictors.  To  be  significant  at  the  99 
percent  level,  a  two  parameter  linear  regression  model  must  explain 
10.3  percent  of  the  variance  and  a  four  parameter  linear  regr^;s1on 
model  must  explain  14.7  percent.  Thus,  both  hypotheses  produce 
statistically  significant  results  showing  that  there  is  another  way 
of  explaining  the  periodicity  observed  In  the  FAI  data. 

A  further  comment  on  the  significance  of  the  pressure  and 
precipitation  analysis  Is  needed.  We  have  seen  the  oscillating 
nature  of  the  pressure  response  at  the  amplitude  modulated  nodal 
frequency.  Standard  statistical  tests  for  regression  look  at  the 
percent  explained  variance  to  calculate  an  F-ratlo  and  from  there  a 
significance  level.  Now,  If  a  station  Is  located  on  the  pivotal 
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axis  of  the  oscillating  anomaly  pattern,  there  will  be  a  minimal 
response  or  amplitude  at  this  location.  Therefore,  the  explained 
variance  will  be  low  and  the  regression  deemed  statistically  not 
significant.  Nevertheless,  the  lack  of  response  at  this  pivotal 
location  Is  physically  significant  and  explainable. 


3.3.3  300  mb  Ug  component. 

In  order  to  enhance  understanding  of  the  synoptic  climatology 
presented  In  the  last  two  sections,  the  response  of  the  300  mb  zonal 
geostrophlc  component  of  the  wind  (Ug)  at  the  modulated  nodal 
frequency  will  be  examined.  Unfortunately  there  were  only  28  years 
of  data  available  (1950-1977)  for  this  analysis  which  allows  only 
1.5  cycles  for  an  18.6  year  oscillation. 

The  Ug  component  was  calculated  using  the  following  equation: 

Ug  -  -(g/f) ah/ay  3.11 

where  f  Is  the  Coriolis  parameter, 

g  Is  the  acceleration  due  to  gravity,  and 
ah/ay  Is  the  north/ south  gradient  of  the  300  mb  heights. 
Since  the  heights  were  available  at  every  5*  of  latitude,  this  was 
the  distance  Interval  used  for  the  height  gradient.  The  value  for  f 
was  calculated  at  the  average  latitude  and  the  values  were  plotted 
at  this  latitude. 

The  percent  variance  map  Is  shown  In  Fig  3.17.  The  nonlinear 
regression  model  for  the  Ug  wind  component  located  at  27.5*N,  75*E 
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Is: 
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14.13Ccos(2«fr1^  ♦  4)  -  l^Dcostevf^t  -  0.25*).  3.12 

Thus  the  phase  (trfth  t«0  In  June  1950)  of  the  forcing  cosine  Is 
alnost  exactly  In  phase  with  the  lunar  nodal  cycle  which  had  a 
Maximum  In  August  1950.  The  correlation  between  this  time  series 
and  the  TMEDLF  response  with  the  susceptibility  term  equal  to  1.04 
Is  0.997.  At  the  same  latitude  the  point  directly  to  the  east  and 
west  had  a  phase  of  1.7*  and  -10.3*,  respectively,  making  them 
nearly  In  phase  with  the  lunar  forcing  also. 

Since  the  early  part  of  the  300  mb  data  record  Is  at  a  time  when 
the  hemisphere  was  near  Its  wannest,  the  amplitude  of  the  response 
Is  quite  low.  Not  until  1970  Is  the  TMEDLF  response  again  achieving 
a  noticeable  fraction  of  Its  maximum  amplitude.  Therefore,  looking 
at  the  anomaly  map  for  1970  (Fig.  3.18)  we  see  a  4  m/s  easterly 
anomaly  between  70*E  and  75*E  at  27.5*N  over  northwestern  India  and 
a  5  m/s  westerly  anomaly  north  of  there  between  65*E  and  70*E.  Both 
anomalies  extend  mainly  eastward  from  the  maximum.  The  last  year  of 
data  Is  1977  which  Is  two  years  before  the  next  peak  In  the  TMEDLF 
response.  Here  we  see  (Fig.  3.18)  the  anomaly  pattern  Is  reversed 
In  sign  with  larger  magnitudes  as  Is  to  be  expected. 

In  order  to  get  a  clear  relationship  between  the  precipitation 
and  300  mb  responses.  It  will  be  necessary  to  extrapolate  the  300  mb 
Ug  component  back  In  time  because  the  11  years  of  overlap 
(1950-1960)  are  at  a  time  when  the  response  Is  minimal.  Since  the 
susceptibility  threshold  for  northern  India  Is  greater  than  one,  the 
300  mb  response  does  not  change  sign  and,  therefore,  all  one  must  do 


Is  look  «t  the  phasing.  The  phasing  Is  such  that  In  1903  a  peak  In 
the  300  Mb  response  would  have  occurred  (as  It  did  In  both  the 
precipitation  and  pressure  responses.)  In  this  year  there  would 
have  been  a  westerly  anomaly.  We  have  already  seen  that  1903  had  a 
dry  precipitation  response  In  the  nodal  cycle.  In  1913  the  Ug 
anomaly  would  have  been  easterly,  and  as  we  have  seen  this  was  a  wet 
year.  Thus,  there  Is  a  westerly  wind  anomaly  In  dry  years  and  vice 
versa.  Since  we  are  dealing  with  a  periodic  forcing  the  validity  of 
this  extrapolation  Is  acceptable;  however.  It  would  be  more 
reassuring  If  a  comparison  could  be  made  In  a  year  covered  by  both 
data  sets.  Therefore,  If  we  look  at  anomalies  from  a  regression  of 
precipitation  that  covers  the  years  1895-1978  we  see  that  In  1978 
when  there  was  a  westerly  anomaly  that  the  precipitation  pattern  was 
dry  (Fig.  3.19)  although  the  anomaly  magnitudes  are  still  relatively 
small. 

3.3.4  Possible  physical  relationship  between  lunar  forcing  and  the 
observed  climatological  data. 

Up  to  this  point,  all  that  has  been  looked  at  Is  the  correlation 
of  the  TMEDLF  response  (based  on  atmospheric  susceptibility  to  lunar 
forcing)  and  the  station  pressure,  precipitation  and  300  mb  Ug 
component.  The  station  pressure  anomaly  patterns  establish  wind 
flow  regimes  that  determine  the  character  of  the  precipitation 
anomalies  and  all  three  of  these  climatological  data  sets  are  In 
phase  with  the  hypothesized  forcing  In  northern  India.  Why? 


To  begin  to  answer  this  question  let  us  apply  the  300  mb 
susceptibility  factor  for  northern  India  to  the  latitudinal 
distribution  of  the  lunar  potential  nodal  cycle  anomaly  as  follows: 

HQgfd)  .  14.1[Q(t)  -  1.04]FQ5(d)  3.13 

where  4  Indicates  a  latitude  dependence.  Then  In  years  when  the 
Ug  component  anomaly  Is  negative  at  the  latitude  of  northern  India 
we  get  a  latitudinal  distribution  for  Hgg  which  has  the  same 
gradient  characteristics  as  seen  In  Figure  2.5.2. 

If  we  were  to  calculate  a  geostrophlc  wind  associated  with  this 
potential  distribution  we  would  see  that  an  easterly  component  Is 
produced.  This  was  the  case  In  1970  when  there  were  stronger  than 
normal  easterlies  south  of  the  Himalayas  and.  In  fact,  since  the 
forcing  and  300  mb  Ug  component  are  In  phase,  every  year  with  this 
potential  anomaly  distribution  has  easterly  anomalies  south  of  the 
Himalayas.  This  Is  a  physically  consistent  relationship. 

However,  to  the  north  of  the  Himalayas  there  are  westerlies  In 
these  years;  and  yet,  as  the  distribution  shows  this  region  Is  also 
Influenced  by  the  same  latitudinal  gradient  of  potential  as  the 
region  south  of  the  Himalayas. 

The  following  possible  explanation  Is  provided.  In  May  and  June 
the  major  synoptic  event  In  the  upper  levels  Is  the  movement  of  the 
200  mb  ridge  axis  to  the  north  from  10*N  to  30*N.  This  movement 
decreases  the  westerlies  In  northern  India,  the  subtropical  jet 
moves  north  of  the  Himalayas  and  easterlies  become  established 
across  northern  India.  In  years  when  this  synoptic  feature  Is 
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subjected  to  an  additional  easterly  component  from  lunar  tidal 
forcing,  the  westerlies  to  the  north  of  the  ridge  weaken  faster  than 
normal  allowing  the  ridge  axis  to  move  northward,  establishing 
easterlies  across  northern  India  earlier.  This  In  turn  allows  the 
subtropical  jet  to  move  north  of  the  Himalayas  earlier  than  normal 
Increasing  the  westerlies  In  that  region.  Thus,  there  Is  an 
easterly  anomaly  to  the  south  of  the  Himalayas  and  a  westerly 
anomaly  to  the  north. 

The  lunar  tidal  forcing  effect  Is  being  seen  here  because  there 
Is  some  singularity  on  which  to  act.  The  ridge  axis  moves  north 
every  year;  however.  In  some  years  Its  northward  movement  Is  impeded 
and  In  other  years  helped  by  the  small  additional  zonal  component 
Induced  by  the  tidal  potential  distribution.  And  as  discussed 
earlier,  the  lunar  forcing  effect  Is  also  contingent  upon  the 
strength  of  the  monsoon  circulation  for  a  given  year.  If  the 
circulation  Is  particularly  strong,  as  It  Is  In  warm  years,  then  the 
effect  Is  not  very  noticeable  at  all. 

The  lunar  forcing  mechanism  acts  on  a  hemispheric  scale  (e.g. 
see  Fig.  2.5)  and  yet  the  observed  response  occurs  on  a  regional 
scale  (l.e. ,  the  northern  half  of  India  In  June).  As  noted  In 
Section  2.3  the  nodal  cycle  (.0537  year'1)  has  only  a  latitude 
dependence.  An  analogous  situation  occurs  with  the  Sun's  thermal 
forcing  which  also  acts  on  a  hemispheric  scaie  with  the  amount  of 
solar  radiation  received  at  the  top  of  the  atmosphere  having  only  a 
latitude  dependence.  Then  because  of  the  difference  In  response  of 
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the  land  and  ocean  to  the  thermal  forcing,  a  strong  temperature 
contrast  develops  producing  the  regional  monsoon  circulation  over 
India.  The  lunar  mechanical  forcing  which  Is  acting  over  the  entire 
Earth  also  produces  a  regional  response  through  an  Interaction  with 
the  monsoon  circulation  and  the  regional  topography*  most  notably 
the  Himalayas. 

Figure  2.5  shows  the  tidal  potential  anomaly  values  for  the 

A  A 

extremes  of  the  nodal  cycle.  The  value  Is  0.0400  ms  at  25*N 
2  ?  • 

and  0.0216  m  s  at  30  N,  giving  a  geostrophlc  zonal  component 
of: 

Ug  -  -Cl/f) ad/ay  -  3.14 

-(0.0216-(0.0400))/(6.73X10“5) (5*) (111000)  *  0.00049  ms*1 
The.,  observed  response  In  the  300  mb  data  Is  almost  10000  times  this 
value  In  1970. 

Thus,  we  are  not  dealing  with  an  equilibrium  tide.  This  Is  not 
surprising,  however,  since  the  ocean  tides  do  not  respond  as  an 
equilibrium  tide  due  to  Irregularities  In  the  ocean  basin  (e.g.,  the 
enormous  tidal  amplitudes  In  the  Bay  of  Fundy).  In  the  region  of 
northern  India  there  Is  the  largest  and  highest  mountain  barrier  In 
the  world  which  must  have  an  Impact  on  the  magnitude  of  the 
atmospheric  lunar  tides.  Nash  (1982)  has  shown  that  the  atmospheric 
pressure  response  to  the  Chandler  pole  tide  Is  also  nonequlllbrlun. 
with  magnitudes  100  to  1000  times  what  Is  predicted  by  equilibrium 
theory.  Currie  (1981)  found  evidence  of  the  nodal  cycle  In 
temperature  and  drought  conditions  and  noted  that  the  atmospheric 
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tide  Is  "strongly  nonequilibrium".  He  also  found  that  the  nodal 
signal  only  appeared  In  regions  which  parallel  the  core  of 
westerlies  In  winter.  Similarly,  as  the  analysis  of  the  monsoon 
circulation  seems  to  Indicate  the  jet  stream  plays  a  major  role  In 
the  Interaction  of  the  lunar  forcing  and  the  atmosphere.  Others 
(Duperfer,  1946  and  Chapman  and  Llndzen,  1970)  have  provided 
evidence  that  the  short  period  lunar  tides  produce  highly 
nonequilibrium  atmospheric  tides  at  higher  altitudes.  Thus,  It  Is 
not  unreasonable  to  assume  that  the  longer  period  tides  behave  In  a 
similar  manner,  producing  the  observed  response  In  the  300  mb  zonal 
geostrophlc  wind  component. 

Since  the  nature  of  this  tide  Is  so  highly  nonequllibrlum,  the 
fairly  straightforward  explanation  put  forth  In  this  section  may 
very  well  be  Invalid.  However,  the  characteristics  of  the  dynamic 
long  period  atmospheric  tides  are  not  In  the  literature;  and  until 
research  In  this  area  Is  performed,  there  Is  no  reference  point  from 
which  to  proceed  with  a  more  dynamically  based  explanation.  Thus, 
for  the  present  time,  the  empirically  deduced  physical  link  between 
forcing  and  response  must  suffice. 

3.4  Response  to  lunar  forcing  at  .263  year-1. 

3.4.1  Introduction. 

After  the  nodal  frequency  Is  removed,  the  frequency  that 
explains  the  most  Interannual  variance  In  the  lunar  potential 
perlodogram  Is  .263  year-1.  The  characteristics  of  the  Moon's 


orbit  that  produce  this  frequency  have  already  been  discussed.  It 
now  remains  to  look  at  the  response  In  the  station  pressure, 
precipitation,  and  300  mb  Ug  component  at  this  frequency. 

As  with  the  nodal  cycle  response  we  will  apply  the  same 
susceptibility  factor  to  the  forcing  frequency.  Further 
justification  for  doing  this  can  be  found  by  looking  at  the  response 
In  the  coefficient  time  series  of  the  first  eigenvector  of  June 
precipitation  (CBB,  Fig.  6a)  over  a  moving  19  year  window.  Each  19 

year  time  series  (e.g.,  1895-1913,  1896-1914 .  1957-1975)  was 

regressed  against  Acos[2»{ .263)t+d].  The  19  year  period  allows  for 
five  cycles  of  the  3.8  year  period.  A  plot  of  the  normalized 
amplitudes  versus  the  last  year  of  the  regression  period  (Fig.  3.20) 
shows  that  this  frequency  was  much  more  Important  In  the  early  and 
latter  parts  of  the  record  than  It  was  In  the  middle.  The 
susceptibility  factor  applied  to  the  forcing  frequency  captures  this 
amplitude  modulation.  By  regressing  the  time  series  against  a 
single  cosine  with  frequency  .263  year"1  one  essentially  gets  an 
average  amplitude  over  the  data  period  (as  suggested  In  Fig.  3.20); 
whereas,  by  regressing  against  the  amplitude  modulated  model  one  can 
capture  the  larger  amplitudes  where  they  occur. 

3.4.2  Station  pressure. 

The  percent  variance  map  (Fig  3.21)  shows  that  the  largest 
response  to  the  nonlinear  model  Is  In  northeastern  India  along  and 
to  the  north  of  the  mean  position  of  the  monsoon  trough.  The  region 
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Figure  3.20.  Normalized  amplitudes  of  the  regression  of  the  first  eigenvector  coe 
(June  precipitation)  against  the  cosine  with  frequency  .263  year-1  f 
lapping  19  year  periods  plotted  against  the  last  year  of  the  period. 


In  southern  India  also  has  an  appreciable  signal.  The 
susceptibility  threshold,  B,  (Fig.  3.22)  In  the  northeast  ranged 
from  .95  to  2.03  Indicating  that  the  F26  forcing  had  some  effect 
at  all  times  and  did  not  change  sign  except  briefly  at  Calcutta 
which  had  the  .95  value.  Note  that  these  values  are  considerably 
higher  than  the  susceptibility  thresholds  for  the  nodal  cycle 
response  (Fig.  3.7).  The  B  term  In  the  southern  peninsula  range., 
from  .35  to  .53. 

The  frequency,  .263  year”*,  produces  a  3.8  year  period; 
thus,  the  anomalies  for  the  first  and  third  years  of  the  data  set 
(1895  and  1897)  provide  a  good  Idea  of  the  pressure  response.  In 
1895  (Fig.  3.23)  there  Is  a  strong  positive  anomaly  extending  all 
the  way  from  the  Bay  of  Bengal  northwestward  along  the  Himalayan 
foothills  to  Lahore  City  (31.6*N,  74.4*E).  There  Is  a  negative 
pressure  anomaly  along  the  west  coast,  which  creates  a  NE/SW 
pressure  gradient.  The  geostrophlcally  balanced  wind  would  blow 
toward  the  northwest.  Such  a  meridional  pressure  gradient  would 
enhance  the  southwesterly  monsoonal  flow;  and  as  we  shall  see  In  the 
next  section,  this  year  Is  wet  through  central  and  northwestern 
India.  Two  years  later  In  1897  the  pattern  has  reversed,  the 
geostrophlc  flow  Is  from  the  northwest  Into  central  India,  and  one 
would  expect  a  dry  year. 

This  pressure  anomaly  pattern  Is  similar  to  the  one  associated 
with  the  response  to  the  lunar  nodal  cycle.  In  both  cases  the  same 
NE/SW  pressure  gradient  Is  established.  Also,  the  maximum  magnitude 


Figure  3.22. 


The  value  of  the  susceptibility  threshold  from  the 
regression  of  the  station  pressure  data  against  the 
kF-c  model  (units  are  value  x  10"2). 
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of  the  anomalies  In  1895  for  both  forcing  frequencies  are  comparable 
(l.e.,  about  1  mb).  One  difference  Is  that  the  western  India 
anomaly  Is  not  as  well  developed  In  the  F26  response  as  It  Is  In 
the  Fq£  response. 

There  Is  another  anomaly  center  located  in  southern  India.  In 
1895  the  negative  anomaly  In  this  region  produces  easterly  flow  to 
Its  north  Inhibiting  the  southwesterly  flow  and  creating  dry 
conditions  along  the  west  coasts  as  Is  confirmed  In  the  next 
section.  In  1897  the  southern  India  high  pressure  anomaly  enhances 
the  southwesterly  flow  to  Its  north  producing  a  positive 
precipitation  anomaly. 

In  1895  the  geostrophlc  wind  anomaly  between  Kota  and  Agra 
(station  numbers  14  and  8,  respectively  In  Fig.  3.4)  Is  2.7  ms"1, 
and  between  Nagpur  and  Allahabad  (numbers  23  and  15)  it  Is  2.2 
ms"1  from  the  southeast  In  a  region  where  the  mean  merlodlonal 
wind  Is  zero  with  standard  deviation  2.1  ms"1  (Crutcher,  1966). 

Because  of  the  amplitude  modulation.  In  1950  (Fig.  3.24)  the 
positive  and  negative  anomaly  centers  are  still  discernible  but 
their  magnitude  Is  greatly  decreased  and  the  low  pressure  anomaly  Is 
not  as  far  to  the  north  against  the  Himalayas. 

To  establish  whether  the  pressure  response  Is  In  phase  with  the 
lunar  forcing  of  similar  frequency,  a  pressure  gradient  time  series 
was  calculated  by  again  taking  the  difference  between  the  anomaly 
average  of  the  eight  northeastern  stations  and  seven  western 
stations.  The  average  susceptibility  threshold  for  these  two  groups 
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of  stations  was  1.20  yielding  a  correlation  between  the  TMEDLF 
response  and  observed  response  of  -0.984. 

Because  the  western  anomaly  Is  not  as  well  developed  In  the 
kF26  response,  the  magnitude  of  the  northeast  minus  west  anomaly 
difference  Is  not  as  great  as  that  for  the  kFQ5  response.  For 
example.  In  1895  (a  year  of  near  maximum  for  both  frequencies),  the 
kFQ5  average  anomaly  difference  Is  1.38  mb  and  that  for  the  kF2g 
model  Is  0.92  mb.  Nevertheless,  these  values  are  comparable  as  one 
might  expect  from  the  T)€DLF  response  perlodogram  of  Fig.  3.2.2. 
Because  the  surface  pressure  responds  to  the  two  lunar  forcing 
frequencies  with  different  susceptibility  thresholds,  the  relative 
magnitude  of  the  response  at  .263  year"1  Is  much  higher  than  the 
perl odogram  of  the  lunar  tidal  potential  (Fig.  3.2.1)  suggests,  for 
the  lunar  tidal  potential  time  series  does  not  account  for  the 
susceptibility  of  the  atmosphere  to  this  lunar  forcing. 

3.4.3  Precipitation. 

The  percent  variance  map  for  precipitation  (Fig.  3.25)  shows  a 
large  response  In  central  and  northern  India.  This  region  lies 
between  the  10  June  and  1  July  average  positions  of  the  monsoon 
front  (see  Fig.  1.1),  suggesting  that  a  large  part  of  the  variance 
might  be  due  to  the  timing  of  the  arrival  of  the  monsoon..  There  Is 
also  an  area  In  southern  India  where  a  significant  amount  of 
variance  Is  associated  with  the  hypothesized  model.  The 
susceptibility  thresholds  are  shown  In  Fig.  3.26. 
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Figure  3.26.  The  value  of  the  susceptibility  threshold  from  the 
regression  of  the  precipitation  data  against  the  kF 
model  (units  are  value  x  10-2). 
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The  anomalies  for  1895,  1897,  1925,  and  1927  are  shown  In  Figure 
3.27.  In  1895  the  moist  southwesterly  flow  produced  positive 
precipitation  anomalies  over  the  central  and  northern  regions  of 
India.  In  the  region  of  the  large  positive  pressure  anomaly  the 
precipitation  anomalies  are  close  to  zero  or  negative.  In  the  south 
the  anomalous  low  pressure  produced  negative  precipitation 
anomalies.  In  1897  the  precipitation  anomalies  are  reversed  as  one 
would  expect  having  already  seen  the  associated  pressure  pattern. 

The  years  1925  and  1927  again  represent  the  two  anomaly  extremes  of 
wet  and  dry,  respectively.  Now,  however,  the  magnitude  of  the 
anomalies  has  been  reduced  due  to  the  amplitude  modulation. 

The  correlation  between  the  average  anomaly  for  the  sixteen 

2 

northern  Indian  stations  with  R  greater  than  or  equal  to  six 
percent  and  the  TMEDLF  response  time  series  Is  -0.983  with  the 
average  susceptibility  threshold  equal  to  1.34  Indicating  that  the 
observed  precipitation  response  Is  again  In  phase  with  the  TMEDLF 
response. 

Even  though  the  percent  variance  Is  low,  the  dissimilar  behavior 
of  Jodhpur  (station  number  11)  Is  of  Interest.  In  1895  when  all 
stations  around  It  have  positive  anomalies,  Jodhpur  has  a  negative 
anomaly  and  vice  versa  in  1897.  Jodhpur  Is  the  "bullseye"  In  Figure 
3.27  In  northern  India.  Bryson  and  Swain  (1981)  reported  that 
rainfall  at  Jodhpur  was  highly  sensitive  to  whether  the  hemisphere 
was  warming  or  cooling  -  that  when  the  hemisphere  was  warming  the 
rainfall  for  July  and  August  was  suppressed.  The  cosine  curve  used 


VOfm’  O**  oWd*  ’•  •  .  ** .  »*  *  V  . '  V  *,*V 


Figure  3.27.  Normalized  anomaly  maps  of  precipitation  for  the  years 
Indicated  from  the  kF?fi  model  (units  are  standard 
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to  approximate  the  hemispheric  termperature  curve  Is  changing  fast 
In  1895  (approximating  warming)  and  rainfall  at  Jodhpur  Is  below 
normal  when  all  other  stations  have  positive  anomalies.  When 
Jodhpur  and  Bikaner  are  regressed  against  a  straight  cosine  with 
frequency  of  .263  year"*,  they  are  both  In  phase.  It  Is  not  until 
the  approximation  of  the  hemispheric  temperature  trend  Is  applied  to 
the  forcing  cosine  that  this  opposing  effect  appears.  In  1925  and 
1927  when  the  rate  of  warming  has  decreased  these  two  stations  are 
In  phase.  Jodhpur's  response  has  changed  sign  because  the 
susceptibility  threshold  Is  0.57  whereas  It  is  0.92  for  Bikaner.  As 
Bryson  and  Swain  (op.  clt.)  showed,  once  the  rate  of  wanning  slowed, 
the  two  stations  had  similar  responses. 

When  discussing  the  kF0S  model  for  precipitation,  studies  by 
Bryson  (1975)  and  Bhalme  and  Mooley  (1981)  were  examined  to  see  how 
well  the  hypothesized  model  fit  the  data  and  as  reported  significant 
variance  was  explained.  Both  studies  involved  precipitation  over 
the  total  monsoon  period.  As  a  result  the  frequency  0.263  year"1 
did  not  appear  as  an  Important  factor  in  explaining  the  Interannual 
variance  In  either  of  these  studies,  for  this  frequency  Is 
associated  with  a  specific  month  and  Is  out  of  phase  with  other 
months.  This  Is  because  .263  year"1  Is  related  to  a  13.7  day 
period  (1/2  the  sidereal  month)  whose  phasing  changes  from  month  to 
month  because  the  length  of  the  calendar  month  Is  not  an  Integral 
number  of  13.7  day  periods  (see  Section  2.4  for  a  further 
explanation).  On  the  other  hand,  the  nodal  period  Is  18.6  years 
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long  so,  for  example,  a  minimum  In  the  cycle  for  one  month's  average 
Is  followed  by  a  minimum  In  the  next  month's  average  (i.e.,  this 
cycle  Is  In  phase  from  month  to  month). 

3.4.4  300  mb  data. 

The  Ug  component  Is  calculated  from  the  300  mb  heights  as  In 
section  3.3.4  using  Equation  3.11.  Regression  of  the  time  series 
for  Ug  at  locations  along  75*E  against  a  cosine  with  frequency  .263 
year"1  yielded  the  following: 

Percent 

Phase  Variance  Amp11tude(ms~l) 

Ug  at  37.5’N  174*  8.3  2.1 

Uj  at  32.5*N  -95*  0.3  0.5 

Ug  at  27.5*11  -3*  10.1  2.0 

After  applying  the  susceptibility  factor  and  using  nonlinear 
regression  we  have: 


Percent  Susceptlbl-  1970 
Phase  Variance  llty  Threshold  Value 

Uq  at  37.5*N  177  9.1  1.3  2.5 
UX  at  32.5  N  -8  7.5  0.8  1.4 
Ug  at  27.5*N  4  10.3  1.8  -2.1 


In  both  regressions  the  Ug  component  at  locations  north  and  south 
of  the  Himalayas  are  out  of  phase.  The  response  of  the  middle 
location  changes  from  what  appears  to  be  a  node  In  the  linear 
regression  to  a  response  that  Is  closely  In  phase  with  either  the 
location  to  the  north  or  to  the  south  of  the  Himalayas  depending  on 
the  sign  of  the  susceptibility  factor  which  changes  because  the 
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susceptibility  threshold  Is  less  than  one.  In  either  case,  however, 
there  Isa  definite  oscillation  north  and  south  of  the  Himalayas. 

Looking  at  the  percent  variance  map  (Fig  3.28)  directly  to  the 
west  of  the  Himalayas  (e.g.,  along  6Q*E)  the  response  Is  weaker  than 
along  the  longitudes  of  the  Himalayas,  and  downstream  an  even 
stronger  response  Is  observed.  For  example,  along  100*E  at  42.5°N 
and  at  27.5 *N  the  nonlinear  model  explains  24.0  percent  and  26.8 
percent  of  the  variance,  respectively.  At  latitude  32.5*N  the 
explained  variance  Is  less  than  it  Is  to  the  north  and  to  the  south 
all  the  way  from  35*E  to  120*E  suggesting  that  a  nodal  line  exists 
In  this  region.  As  we  saw  In  Chapter  2  there  Is  a  node  at  35.26*N 
In  the  lunar  potential.  Like  the  response  along  75*E,  the  response 
downstream  of  the  Himalayas  Is  out  of  phase  north  and  south  of  this 
nodal  line  as  seen  In  the  anomaly  pattern  for  1971  (Fig  3.29). 

These  results,  along  with  similar  ones  for  the  kFQ5  model,  suggest 
that  a  dynamical  Interaction  Involving  a  forcing  at  these  two 
frequencies,  the  Himalayas,  and  the  subtropical  jet  creates  a 
response  In  the  vicinity  of  the  Tibetan  Plateau  which  then 
progaga tes  downstream. 

That  the  observed  Ug  response  at  27.5*N,  75*E  Is,  In  fact.  In 
phase  with  the  kF25  TMEDLF  response  Is  seen  from  their  -.978 
correlation.  Comparing  the  phasing  of  the  precipitation  and  Ug 
response  Indicates  that  In  dry  years  there  Is  a  westerly  Ug 
anomaly  over  northern  India  and  an  easterly  anomaly  north  of  the 
Himalayas.  This  Is  consistent  with  the  results  for  the  kF0g  model. 


Figure  3.28.  The  percent  of  the  geostrophic  zonal  wind  (U  ) 
variance  explained  by  the  kF«g  model  (units  9 
are  percent).  Level  is  300  millibars. 


3.4.5  Possible  physical  relationship  between  the  hypothesized  lunar 

forcing  and  the  observed  climatological  data. 

As  with  the  kFQg  model.  If  we  apply  the  susceptibility  factor 
from  the  300  mb  Ug  component  model  In  northern  India  at  27. 5* N, 

75*E  to  the  lunar  tidal  potential  as  follows: 

H26(rf)  »  -1.9[Q(t)  -  1.81]F26(tf)  3.15 

then  In  years  when  the  Ug  component  anomaly  Is  negative  (as  In 
1971,  Fig.  3.29)  at  this  location,  the  latitudinal  distribution  of 
H26  Is  the  same  as  It  was  for  HQ5  In  section  3.3.4.  This 
distribution  would  then  produce  a  negative  or  easterly  U 

5# 

component  If  we  consider  geostrophlc  equilibrium.  Thus,  the 
synoptic  relationships  for  the  kFQ5  response  discussed  earlier 
apply  equally  well  for  the  kF26  response,  for  both  lunar  forcing 
frequencies  produce  similar  phasings  In  all  three  climatic 
variables. 

3.5  Summary. 

As  we  have  seen,  all  three  climatic  variables  (l.e.,  pressure, 
precipitation,  and  300  mb  Ug)  are  highly  correlated  with  the  lunar 
forcing  response  models,  kFQg  and  kF26,  which  means  they  are 
highly  correlated  with  each  other.  Table  3.4  provides  a  summary  of 
these  statistics.  The  average  susceptibility  threshold  values 
summarized  In  this  table  Indicate  that  the  values  used  to  calculate 
the  perlodogram  of  Fig.  3.2.2  were  quite  reasonable. 

As  we  have  seen,  a  low  pressure  anomaly  along  the  foot  of  the 
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Table  3.4 

Summary  of  statistics  for  the  kFgg  and  kFgg  models 
for  the  three  climatic  variables 


Average 
Susceptibi¬ 
lity  threshold 


Average 
Susceptibi¬ 
lity  threshold 


Himalayas  plus  a  high  pressure  anomaly  along  the  west  coast  leads  to 
deficient  rainfall  over  most  of  India.  Ramaswamy  (1958)  examined 
the  mean  zonal  wind  component  at  New  Delhi  and  Calcutta  for  both  a 
weak  and  active  monsoon  period  and  found  that  at  both  stations  the 
dry  period  was  characterl zed  not  only  by  a  northward  shift  of  the 
low  pressure  trough  to  a  position  along  the  Himalayas  but  also  by  a 
westerly  anomaly  at  these  stations  up  to  very  high  levels.  More 
recently  Raman  and  Rao  (1981)  have  shown  that  establishment  of  an 
east  Asia  blocking  ridge  along  100*E,  extending  from  35*N  to  70*N, 
anchors  a  westerly  trough  around  70*E  to  75*E.  With  this  trough 
extending  Into  more  southerly  latitudes  (30*N)  than  normal  the  lower 
tropospheric  monsoon  trough  Is  then  displaced  northward  into  the 
Himalayas.  This  synoptic  situation  would  produce  a  westerly  anomaly 
at  300  mb  In  the  region  where  we  see  It  In  Figure  3.18  for  the  year 


Between  April  and  July  along  75*E  the  subtropical  jet  moves  from 
25*N  to  45*N  while  the  200  mb  ridge  axis  moves  from  10*N  to  30*N. 

Yin  (1949),  Wright  (1967)  and  many  others  have  found  that  the  timing 
of  the  onset  of  the  monsoon  Is  directly  related  to  the  movement  of 
the  subtropical  jet  from  Its  winter  postion  south  of  the  Tibetan 
Plateau  to  Its  summer  position  to  the  north  of  the  plateau.  Wright 
(op.  clt.)  derived  several  relationships  between  the  200  mb  zonal 
wind  and  the  onset  of  the  monsoon.  He  found  that  the  first  main 
burst  of  the  monsoon  extended  up  to  at  least  13*N  and  occurred 
within  "a  fortnight  of  the  date  on  which  the  200  mb  westerly  wind 


component  at  Bombay  falls  permanently  below  5  kt  For  the 

period  1891-1950  the  average  date  for  the  monsoon  to  arrive  at  13*N 
was  3  June  with  a  standard  deviation  of  6.3  days.  The  second  main 
burst  of  the  monsoon  extended  up  to  22*N  and  occurred  "about  a  month 
after  the  start  of  the  main  decrease  In  the  200  mb  westerly  wind 
component  at  New  Delhi  ...".  The  arrival  date  at  22°N  was  20  June 
with  a  range  of  22  days.  The  observed  climatic  response  to  both 
hypothesized  lunar  forcing  frequencies  has  shown  that  when  a 
westerly  anomaly  at  300  mb  lies  across  northern  India,  a  low 
pressure  anomaly  lies  along  the  Himalayas  and  rainfall  Is 
deficient.  This  Is  consistent  with  Wright's  analysis  since  late 
arrival  of  the  monsoon,  and  hence  below  normal  rainfall,  was 
associated  with  stronger  than  normal  westerlies.  It  Is  Important  to 
note  that  the  precipitation  appears  to  be  related  more  to  the 
strength  of  the  westerly  anomaly  over  northern  India  than  to  the 
actual  position  of  the  subtropical  jet.  For  the  27  years  of  300  mb 
data,  the  jet  was  never  located  between  25*N  and  30*N  for  the  month 
of  June  and  yet  It  Is  the  zonal  component  of  the  geostrophlc  wind 
anomaly  In  this  region  that  Is  so  highly  correlated  to  the 
precipitation  at  these  two  lunar  frequencies. 

From  the  analysis  In  this  chapter,  then,  a  consistent 
climatology  emerges  that  Is  similar  for  both  lunar  forcing  frequency 
models: 

a.  There  Is  an  external  boundary  condition,  lunar  tidal 
forcing,  that  Is  always  acting  upon  the  earth-atmosphere-ocean 
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system.  This  Is  an  ultimate  cause  of  the  climatic  state  which 
establishes  an  array  of  weather  patterns  associated  with  the  monsoon 
circulation. 

b.  South  of  the  Himalayas*  the  300  mb  geostrophlc  zonal 
component  of  the  wind  responds  to  the  latitudinal  distribution  of 
the  lunar  potential  anomalies  resulting  from  the  two  frequencies, 
.0537  year”1  and  .263  year”1.  These  two  frequencies  account  for 
almost  100  percent  of  the  lunar  tidal  potential  Interannual 
variability  In  June.  At  one  extreme  of  each  cycle  there  Is  a 
westerly  geostrophlc  wind  anomaly  and  at  the  other  extreme  there  Is 
an  easterly  anomaly. 

c.  In  years  when  there  Is  a  westerly  anomaly  at  27.5*N  over 
northern  India,  the  normal  decrease  In  westerlies  at  this  latitude 
Is  slowed.  As  other  research  has  shown,  the  timing  of  the  onset  of 
the  monsoon  Is  directly  related  to  the  timing  of  the  decrease  In  the 
westerlies  over  northern  India.  Thus,  a  westerly  anomaly  In  this 
region  results  In  a  delay  of  the  onset. 

d.  This  delay  In  the  monsoon  onset  Is  related  to  the  surface 
pressure  anomaly  pattern  that  occurs  at  the  same  time  as  the  300  mb 
westerly  anomaly.  This  surface  pressure  anomaly  pattern  Is  one  of 
low  pressure  along  the  base  of  the  Himalayas  and  high  pressure  In 
western  India  producing  a  NE/SW  pressure  gradient  that  supports 
anomalous  northwesterly  flow  In  direct  opposition  to  the  movement  of 
the  monsoon  front  toward  the  northwest.  The  high  pressure  anomaly 
along  the  west  coast  also  decreases  the  southwesterly  monsoonal  flow 
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off  of  the  Arabian  Sea.  With  the  normal  low  level  monsoon 
circulation  disrupted  In  this  manner,  the  onset  Is  delayed. 

e.  In  northern  India  most  of  the  June  rainfall  occurs  after  the 
monsoon  has  arrived.  The  normal  onset  date  Is  late  June. 

Therefore,  aqy  delay  In  this  date  will  have  a  large  Impact  on  the 
percent  of  normal  rainfall  received  In  June. 

f.  As  CBB  showed,  the  precipitation  response  at  the  two  lunar 
tidal  frequencies,  .263  year**  and  .0482  year'1  (the  linear 
signature  of  the  nodal  cycle  response).  Is  observed  to  be  most 
Important  In  explaining  the  Interannual  variability  of  precipitation 
In  northern  India  In  June.  Thus,  It  would  appear  that  lunar  tidal 
forcing  plays  a  significant  role  In  the  Interannual  variability  of 
June  precipitation  In  northern  India. 

g.  The  Northern  Hemisphere  surface  temperature  Is  a  proximal 
cause  of  the  climate,  for  It  is  related  to  the  amount  of  solar 
radiation  reaching  the  Earth's  surface  which  Is  controlled  by  such 
ultimate  causes  as  volcanic  activity  and  the  Earth's  orbital 
parameters.  The  Northern  Hemisphere  surface  temperature  determines 
to  a  large  extent  the  Intensity  of  the  monsoon  circulation.  As  a 
result,  this  proximal  cause  modulates  the  effects  of  the  lunar  tidal 
forcing  described  above.  When  the  hemisphere  Is  warm,  the  monsoon 
circulation  Is  strong  and  the  atmosphere  Is  less  susceptible  to  an 
external  forcing,  thus,  diminishing  the  effect  of  the  lunar  tide.  A 
nonlinear  regression  model  was  developed  based  on  this  Idea  of 
atmospheric  susceptibility.  As  noted  at  the  time  the  susceptibility 


factor  used  was  an  approximation.  Nevertheless,  applying  this 
factor  to  the  hypothesized  forcing  frequencies  virtually  "unlocked" 
the  nodal  cycle  In  the  data  and  made  the  amplitudes  of  the  response 
at  the  beginning  of  the  record  for  both  lunar  frequencies  at  least 
two  to  three  times  larger  than  those  attained  In  a  straight  linear 
regression. 


CHAPTER  4 


Conclusions 

1.  This  research  has  presented  definitive  evidence  that  there 
Is  an  Important  periodic  component  In  the  Interannual  variability  of 
pressure,  precipitation,  and  upper  air  winds  for  the  monsoon  onset 
month  of  June  and  that  the  two  dominant  frequencies  found  In  the 
precipitation  records  of  northern  India  stations  are  the  same  two 
frequencies  which  explain  almost  100  percent  of  the  Interannual 
variance  of  the  hypothesized  forcing  mechanism  -  monthly  mean  lunar 
tidal  potential.  Not  only  are  the  response  and  forcing  frequencies 
the  same  but  they  are  also  In  phase  which  strongly  suggests  that 
lunar  tides  In  the  atmosphere  do.  In  fact,  produce  an  element  of 
climatic  variability. 

2.  Past  research  has  typically  examined  monsoon  variability 
using  a  stratified  climatology  approach,  l.e.,  by  studying  certain 
variables  during  drought,  flood,  and/or  normal  years.  However,  by 
examining  the  periodic  component  of  climatic  variability  we  can 
Include  drought,  flood  and  normal  years  Into  a  single  model  based  on 
a  periodic  forcing  mechanism.  Both  lunar  forcing  frequencies, 
examined  In  this  research,  are  associated  with  consistent  climatic 
responses  in  surface  pressure,  precipitation  and  the  300  mb 
geostrophlc  zonal  wind  component.  As  a  result,  the  work  of  many 
monsoon  researchers  could  be  synthesized  Into  a  coherent  synoptic 


climatology  defined  by  a  nonlinear  periodic  model  of  climatic 
varl  ability. 

3.  A  trend  In  the  amplitude  of  the  periodic  response  was  found 
and  related  to  the  Northern  Hemisphere  surface  temperature.  Thus* 
In  this  case  the  trend  and  periodic  components  of  climatic 
variability  are  Inseparable.  This  dependence  of  the  periodic 
component  on  the  hemispheric  temperature  Indicates  that  the  period 
of  time  covered  by  a  data  set  can  determine  whether  or  not  an 
hypothesized  periodic  forcing  will  produce  an  observable  slgn&l  In 
the  data. 

4.  The  synoptic  climatology  associated  with  the  responses  of 
surface  pressure,  precipitation  and  the  300  mb  Ug  component  were 
similar  for  both  lunar  forcing  frequencies.  However,  the 
characteristics  of  the  susceptibility  factors  at  a  given  location 
for  each  frequency  can  be  quite  different.  For  example,  the 
susceptibility  threshold  Is  generally  higher  for  the  kF2g  response 
than  that  for  kFgj.  As  we  saw,  this  resulted  In  the  magnitude  of 
the  k F26  response  being  comparable  to  that  for  the  kF^ 
response.  A  comparison  of  Eqns.  3.13  and  3.15  Indicates  another 
difference  In  that  the  signs  on  the  susceptibility  factor 
coefficients  are  opposite;  and  since  In  both  cases  the  [Q(t)-B] 
factor  Is  always  less  than  zero,  the  sign  on  the 


forcing  factor  Is  negative  for  the  kFQg  response  and  positive  for 
the  kFgg  response.  Thus,  opposite  latitudinal  distributions  of 
the  lunar  tidal  potential  produce  similar  300  mb  Ug  anomalies.  In 
a  monthly  average  we  can  assume  geostrophlc  equilibrium  with  the 
corlolls  force  (C)  balancing  the  pressure  gradient  force  (P). 
However,  with  the  hypothesis  of  lunar  tidal  forcing,  there  Is  an 
additional  force  (L)  to  consider.  With  geostrophlc  equilibrium,  we 
have. 

Cap.  4.1 

Considering  the  sign  on  the  forcing  for  the  kFg5  response  we  have, 

C  ■  P  -  L  >  P  ■  C  +  L  ;  4.2 

and  for  the  kF26  response  we  have, 

C-P  +  L>P  +  L-  C.  4.3 

The  reason  for  this  difference  In  sign  Is  not  known  at  the  present 
time  but  It  may  be  due  to  the  different  time  scales  associated  with 
each  frequency.  f26  arises  from  a  13.7  day  period  which  goes 
through  two  complete  cycles  during  a  calendar  month  with  a  2.7  day 
residual  In  the  month  of  June.  On  the  other  hand,  F05  arises  from 
the  nearly  6800  day  nodal  cycle  and,  therefore,  acts  In  the  same 
sense  throughout  the  entire  calendar  month.  Such  a  difference  In 
time  scale  certainly  has  Implications,  for  example.  In  the 
Interaction  between  the  atmosphere  and  the  oceans. 

5.  The  Impetus  for  this  research  was  to  Improve  long  range 
climate  forecasting  capabilities.  This  research  has  provided  an 


analysis  of  one  segment  of  the  past  climatic  variability  during  the 
monsoon  onset  month  of  June  and  found  that  there  Is  an  important 
periodic  component  to  this  variability  that  Is  In  phase  with  lunar 
tidal  forcing.  Therefore,  there  Is  an  inherent  predictive 
capability  In  the  nonlinear  model  used,  for  the  Moon's  position  can 
be  very  accurately  predicted.  The  other  factor  In  the  model. 
Northern  Hemisphere  surface  temperature,  must  also  be  forecast. 

These  forecasts  are  currently  available  using  an  atmospheric  optical 
depth  model.  The  susceptibility  factor  must  then  be  modified  to  use 
the  temperature  forecast  Instead  of  the  cosine  approximation  used  in 
this  research. 

6.  The  nature  of  the  nonlinear  kF0g  response  shifts  the 
spectral  peak  In  the  forcing  at  .054  year"1  to  a  lower  frequency 
around  .048  year-1.  As  we  have  seen,  this  Is  very  close  to  the 
sunspot  frequency  and  both  hypothesized  forcings  happen  to  be 
closely  In  phase  during  the  period  of  time  covered  by  the  data  sets 
used  In  this  research.  The  lunar  tidal  forcing  hypothesis, 
therefore,  sheds  some  doubt  on  the  sunspot  hypothesis  at  least  to 
the  extent  that  this  research  provides  another  explanation  for  the 
periodicity;  and  certainly,  at  this  point  In  time,  there  is  a 
stronger  physical  relationship  between  atmospheric  tides  and  the 
observed  climatic  variability  than  there  Is  for  the  sunspot 
hypothesis.  This  finding  also  points  out  the  limitations  of  a 
linear  analysis  like  the  Fast  Fourier  Transform  or  perlodogram,  for 


the  true  forcing  frequency  will  not  necessarily  show  up  If  It  Is 
being  modulated  by  some  susceptibility  factor. 

7.  The  averaging  period  of  the  data  used  In  one's  research  will 
determine  what  lunar  signals  are  observable.  For  example,  data 
averaged  over  a  27.3  day  period  would  not  show  a  signal  at  .263 
year"*.  Thus,  the  meteorologist's  standard  method  of  averaging 
over  a  30  or  31  day  period  allowed  this  frequency  to  show  up  In  the 
data.  Similarly,  an  average  of  June,  July,  August,  and  September 
data  would  wash  out  the  .263  year"*  signal  since  this  frequency 

has  a  different  phase  In  each  month.  On  the  other  hand  because  the 
nodal  cycle  Is  such  a  long  period  phenomenon  any  averaging  period  at 
least  up  to  a  year  would  allow  It  to  appear. 

8.  External  forcing  on  the  atmosphere  seems  to  be  observable  on 
the  Interannual  time  scale  only  In  regions  where  a  synoptic 
singularity  exists  (e.g.,  the  arctic  front  In  December  In  the  upper 
midwest  or  the  onset  of  the  monsoon  In  June  In  northern  India)  and 
then  only  If  the  atmosphere  Is  susceptible  to  that  forcing. 
Therefore,  these  sensitive  regions  need  to  be  Identified  so  that  the 
predictive  capability  of  periodic  external  forcing  can  be  explored 
further. 


REFERENCES 


Adderl oy,  E.E.  and  E.G.  Bowen,  1962:  Lunar  component  In 
precipitation  data.  Science.  Vol  137,  749-750. 

American  Ephemerls  and  Nautical  Almanac  (1895-1960),  United  States 
Government  Printing  Office,  Washington,  D.C.,  1980 

Bhalme,  H.N.,  1972:  Trends  and  quasi -biennial  oscillation  In  t*»e 
series  of  cyclonic  disturbances  over  the  Indian  region.  1.  iG, 

Vol  23,  355-358.  "  ~ 

_  and  D.A.  Hooley,  1980:  Large-scale  droughts/f 1 oods  and  monsoon 

circulation.  MWR,  Vol  108,  1197-1211. 

_ and _ ,  1981:  Cyclic  fluctuations  In  the  flood  area  and 

relationship  with  the  double  (Hale)  sunspot  cycle.  JAM,  Vol  20, 
1041-1048. 


_ p _ ,  and  S.K.  Jadhav,  1983:  Fluctuations  In  the 

drought/fl ood  area  over  India  and  relationships  with  the 
Southern  Oscillation.  MWR,  Vol  111,  86-94. 

Bhargava,  B.N.  and  R.K.  Bansal,  1969:  A  quasi -biennial  oscillation 
in  precipitation  at  some  Indian  stations.  IJMG,  Vol  20,  127-128. 

Bradley*  D.A.,  1964:  Tidal  components  In  hurricane  development. 
Nature,  Vol  204,  136-138. 

Brier,  G.W.,  1968:  Long-range  predictions  of  the  zonal  westerlies 
and  some  problems  In  data  analysis.  Rev.  Geophys.,  Vol  6, 
525-551. 

_ _ and  D.A.  Bradley,  1964:  The  lunar  synodical  period  and 

precipitation  In  the  United  States.  JAS,  Vol  21,  386-395. 

Brlnkmann,  W.A.R.,  1979:  Associations  between  temperature  trends. 
Annals  of  the  Association  of  American  Geographers,  Vol  69, 
75G3BT. - - 

Brouwer,  D.  and  G.M.  Clemence,  1961:  Methods  of  Celestial  Mechanics. 
Academic  Press,  New  York,  598  pp. 

Brownlee,  K.  A.,  1965:  Statistical  Theory  and  Methodology.  John 
Wiley  a  Sons,  Inc.,  toew  York,  590  pp. 

Bryson,  R.A. ,  1948:  On  a  lunar  bl -f ortnl ghtly  tide  In  the 
atmosphere.  Trans.  Am.  Geophys.  Union.  Vol  29,  473-475. 


Bryson,  R.A. ,  1975:  The  lessons  of  climatic  history.  Environmental 
Conservation.  Vol  2,  163-170. 

_  and  T.B.  Starr,  1977:  Chandler  tides  in  the  atmosphere.  JAS, 

7ol  34,  1975-1986.  - 

_ _ and  B.M.  Goodman,  1980:  Volcanic  activity  and  climatic  change, 

Science,  Vol  207,  1041-1044. 

and  A.M.  Swain,  1981:  Holocene  variations  of  monsoon  rainfall 
Tn  Rajasthan.  Quaternary  Research,  Vol  16,  135-145. 

Campbell,  W.  H.,  J.B.  Blechman,  and  R.A.  Bryson,  1983:  Long-period 
tidal  forcing  of  Indian  monsoon  rainfall:  An  hypothesis.  JCAM, 
Vol  22,  287-296. 

Chapman,  S.  and  R.L.  Lindzen,  1970:  Atmospheric  Tides.  Gordon  and 
Breach,  Hew  York,  200  pp. 

Crutcher,  H.  L.,  1966:  Components  of  the  1000  mb.  winds  (or  surface 
winds)  of  the  Northern  Hemisphere.  NAYAIR  50-lc-51,  Office  of 
Chief  of  Naval  Operations,  77  pp. 

Currie,  R.C.,  1981:  Evidence  for  18.6  year  MN  signal  In 

temperature  and  drought  conditions  In  North  America  since  AD 
1800.  J.,  Geophys.  Res.,  86,  11055-11064. 

Das,  P.K.,  1968:  The  Monsoon.  National  Book  Trust,  New  Delhi,  162  pp 

Doodson,  A.T.,  1921:  The  harmonic  development  of  the  tide  generating 
potential.  Proc.  Roy.  Soc.,  Ser.  A,  Vol  100,  305-329. 

Duperler,  A.,  1946:  A  lunar  effect  on  cosmic  rays.  Nature,  Vol  157, 
296-299. 


Gill,  A.  E.,  1980:  Some  simple  solutions  for  heat-induced  tropical 
circulation.  QJRMS,  Vol  106,  447-462. 

Godbole,  R.V.,  1973:  Nunerical  simulation  of  the  Indian  summer 
monsoon.  IJMG.  Vol  24,  1-14. 

Goldstlne,  H.  H.,  1973:  New  and  Full  Moons  1001  B.C.  to  A.D,  1651. 
American  Philosophical  Society,  Philadelphia,  221  pp. 

Groveman,  B.S.  and  H.E.  Landsberg,  1979:  Reconstruction  of  Northern 
Hemisphere  temperature:  1579-1880.  Publ.  No.  79-181,  June,  1979, 
Meteorology  Program,  Unlv.  of  Maryland,  College  Park,  MD  20742, 
46  pp. 


Hahn,  D.G.  and  S.  Manabe,  1975:  The  role  of  mountains  In  the  south 
Aslan  monsoon  circulation.  JAS.  Vol  32,  1515-1541. 

_  and  J.  Shukla,  1976:  An  apparent  relationship  between  Eurasian 

snow  cover  and  Indian  monsoon  rainfall.  JAS,  Vol  33,  2461-2462. 

Jagannathan,  P.  and  B.  Parthasarathy,  1973:  Trends  and  periodicities 
of  rainfall  over  India.  MWR,  Vol  101,  371-375. 

Jones,  P.  0.,  T.M.L.  Wlgley,  and  P.M.  Kelley,  1982:  Variations  In 
surface  air  temperatures:  Part  1.  Northern  Hemisphere, 

1881-1980.  MWR,  Vol  110,  59-70. 

Koteswaram,  P.  and  S.M.A.  Alvl,  1959:  Secular  trends  and 

periodicities  In  rainfall  at  west  coast  stations  In  India. 
Current  Science.  Vol  38,  229-231. 

Kutzbach,  J.E.,  1970:  Large-scale  features  of  monthly  mean  Northern 
Hemisphere  anomaly  maps  of  sea  level  pressure.  MWR,  Vol  98, 
708-716. 

and  B.  L.  Otto-Bllesner,  1982:  The  sensitivity  of  the 
Afrlcan-Aslan  monsoonal  climate  to  orbital  parameter  changes  for 
9000  Years  B.P.  In  a  low- resolution  general  circulation  model. 
JAS,  Vol  39,  1177-1188. 

Krl shnamurtl ,  T.N.,  1971:  Observational  study  of  the  tropical  upper 
tropospheric  motion  field  during  the  Northern  Hemisphere  summer. 
JAM,  Vol  10,  1066-1096. 

_ ,  S.M.  Daggupaty,  J.  Fein,  M.  Kanaraltsu,  and  J.D.  Lee,  1973: 

“Tibetan  high  and  upper  tropospheric  tropical  circulations  during 
northern  summer.  BAMS,  Vol  54,  1234-1249. 

and  H.N.  Bhalrae,  1976:  Oscillations  of  a  monsoon  system.  Part 
I.  Observational  aspects.  JAS,  Vol  33,  1937-1954. 

Maksimov,  I.V.,  1958:  Nutatlonal  phenomena  In  the  high  latitude 
atmosphere  andthelr  role  In  the  formation  of  climate.  Ottawa, 
Problems  of  the  North,  No  1  (1960),  Hat.  Res.  Council  Canada. 

nn-izr. - 

.  V.P.  Karklln,  E.I.  Sarukhanlan,  and  N.P.  Smirnov,  1967: 
Nutatlonal  migration  of  the  Icelandic  low.  Dokl.  Akad.  Nauk. 

SSSR.  Vol  177,  88-91.  (English  trans.  Dokl.  Acad.  sci.  U55R7 
tirtn  Scl.  Section) 


Matsuno,  T.,  1966:  Quasl-geostrophlc  motions  In  the  equatorial  area. 
JMSJ,  Vol  44,  25-43. 

Nash,  E.  R.,  1982:  A  critical  comparison  of  two  theories  of  the 
atmospheric  Chandler  “pole  tide".  Masters  thesis.  University  of 
Wisconsin,  28  pp. 

Overland,  J.E.,  and  R.W.  Prelsendorfer,  1982:  A  significance  test 
for  principal  components  applied  to  a  cyclone  climatology.  Mon. 
Wea.  Rev..  110,  1-4. 

Parthasarathy,  B.  and  O.A.  Mooley,  1978:  Some  features  of  a  long 
homogeneous  series  of  Indian  summer  rainfall.  MWR,  Vol  106, 
771-781. 

Pearson,  D.  C.,  1982:  Periodic  Interannual  variations  of  mldwestern 
United  States  temperatures  In  December.  Masters  Thesis, 
University  of  Wisconsin,  76  pp. 

Raghavan,  K. ,  1973:  Break-monsoon  over  India.  MWR ,  Vol  101,  33-43. 

Raman,  C.  R.  V.  and  Y.  P.  Rao,  1981:  Blocking  highs  over  Asia  and 
monsoon  droughts  over  India.  Nature,  Vol  289,  271-273. 

Ramaswamy*  C.,  1958:  A  preliminary  study  of  the  behaviour  of  the 
Indian  southwest  monsoon  In  relation  to  the  westerly  jet- stream. 
Geophyslca,  Vol  6,  455-477. 

Rao,  Y.P.,  1976:  Southwest  Monsoon.  India  Met.  Dept.,  Poona,  367  pp. 

Shukla,  J.,  1975:  Effect  of  Arabian  sea-surface  temperature  anomaly 
on  Indian  summer  monsoon:  A  numerical  experiment  with  the  GFDL 
model.  JAS,  Vol  32,  503-511. 

_  and  B.M.  Mlsra,  1977:  Relationships  between  sea  surface 

temperature  and  wind  speed  over  the  central  Arabian  sea,  and 
monsoon  rainfall  over  India.  MWR,  Vol  105,  998-1002. 

Trenberth,  K.E.  and  D.A.  Paollno,  Jr.,  1980:  The  Northern  Hemslphere 
sea-level  pressure  data  set:  Trends,  errors  and  discontinuities. 
MWR,  Vol  108,  855-872. 

van  Loon,  H.  and  J.  Williams,  1976:  The  connection  between  trends  of 
mean  temperature  and  circulation  at  the  surface:  Part  I.  Winter. 
MWR.  Vol  104,  365-380. 

Vlsvanathan,  T.R.,  1966a:  Formation  of  depressions  In  the  Indian 
seas  and  lunar  phase.  Nature,  Vol  210,  406-407. 


,  1966b:  Heavy  rainfall  distribution  In  relation  to  the  phase 
of  the  moon.  IJMG,  Vol  17,  369-372. 

Voigt,  H.  H.,  1974:  Outline  of  Astronomy.  Vol  I..  Moordhoff 
International  Publ.,  Leydon,  The  Netherlands,  272  pp. 

Walker,  G.T. ,  1924:  Correlation  In  seasonal  variation  of  weather  IX 
Mew.  India  Meteor.  Dept.,  Vol  24,  275-332. 

Weare,  B.C.,  1979:  A  statistical  stucjy  of  the  relationships  between 
ocean  surface  temperatures  and  the  Indian  monsoon.  JAS,  Vol  36, 
2279-2291. 

Webster,  P.J.  and  L.C.  Chou,  1980:  Seasonal  structures  of  a  simple 
monsoon  system.  JAS,  Vol  37,  354-367. 

Wood,  Fergus  J.,  1976:  The  Strategic  Role  of  Perlgean  Spring  Tides. 
U.S.  Department  of  Comaerce,  National  bceanlc  and  Atmospheric 
Administration,  Washington,  D.C.,  538  pp. 

Wright,  P.B.,  1967:  Changes  In  200  mb  circulation  pattern-  related 
to  the  development  of  the  Indian  south-west  monsoon.  Met.  Mag., 
96,  302-315.  ”  ~ 

Yin,  M.T. ,  1949:  A  synoptlc-aerologlc  stuty  of  the  onset  of  the 

summer  monsoon  over  India  and  Burma.  Journal  of  Meteorology.  Vol 
6,  393-400. 


END 

FILMED 

1-84 

\ 

DTIC 


